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|  Laser-damage  thresholds  are  tabulated. 


StCu"irv  or  1 


• -i,r  t  ■^*^-'***^^»**.< 


r 


TABLE  OF  CONTENTS 

page 


Abstract  5 

A.  Introduction  7 

B.  Scattering  by  Macroscopic  Imperfections  13 

C.  Defect  Chemistry  of  Alkali  Halides  and  Their  Optical  Properties  31 

D.  Light  Scattering  by  Point  Defects  in  Insulating  Crystals  70 

E.  Multiple  Scattering  86 

F.  Tabulation  of  Infrared-Absorption  Coefficients  101 

G.  V  Curves  of  Absorption  and  Scattering  170 

H.  Extrinsic  Absorption  179 

I.  Tabulation  of  Laser-Breakdown  Thresholds  of  Solids  194 

J.  Optimum  Frequency  Range  201 

K.  Choice  of  Glass-,  Polycrystalline-,  or  Single-Crystal-Fibers  206 

L.  Publications  212 

M.  Summary  and  Conclusions  213 


I 

I 


ABSTRACT 


Extrinsic  scattering  coefficients  (with  values  as  great  as  1  to  100 

dB/km  for  any  one  of  seven  scattering  processes)  and  extrinsic  absorption 

coefficients  (with  impurity  concentrations  of  either  isolated  impurity 

ions  or  molecules  or  macroscopic  inclusions  as  low  as  10"  3  to  10"4  parts 

per  billion  giving  ggb  =  10  dB/km)  must  be  greatly  reduced  in  order 

to  attain  an  extinction  coefficient  B  =  10‘2  dB/km  =  2  x  lo“8  cm"1, 

which  is  extremely  low,  but  is  a  fairly  modest  goal  for  low-loss  fibers. 

There  are  only  a  few  candidate  glasses,  and  current  fabrication  techniques 

of  crystalline  fibers  give  extinction  as  great  as  five  orders  of  magnitude 
_2 

greater  than  10  dB/km.  Attaining,  maintaining,  and  even  measuring  the 
high  purity  in  both  crystals  and  glasses  is  expected  to  be  difficult. 
Impurity-absorption  spectra  and  peak-absorption  wavelengths  of  both 
infrared  and  visible/ultraviolet  absorption  are  tabulated  for  use  in 
identifying  troublesome  Impurities.  Experimental  absorption 
coefficients  in  the  infrared  region  are  tabulated  to  establish  the 
state  of  the  art.  Single-crystal  alkali  halides  may  be  rendered  useless 
as  low-loss  fiber  materials  by  surface  imperfections  (steps)  and  fracture, 
both  induced  by  bending.  Plots  of  intrinsic-scattering-  and  absorption- 
coefficients,  the  so-called  V  curves,  which  are  compiled  here,  are 
useful  in  establishing  the  wavelength  region  in  which  intrinsic  loss 
is  negligible.  But  it  is  not  necessary  to  operate  at  the  minimum 
of  the  V  curve  because  extrinsic  scattering  and  absorption  are  expected 
to  limit  the  loss.  Values  of  the  zero-dispersion-crossover  wavelength 
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2  2 

XQd  (at  which  d  n/dX  =  0)  are  tabulated  because  informati on-carrying 

fibers  must  be  operated  near  x^g*  The  choice  of  the  optimum  operating 

wavelength  involves  compromises  in  view  of  the  conflicts  in  the  above 

requirements.  The  Issues  involved  in  the  choice  between  si ngl e-crystal - 

and  polycrystalline-fibers  are  identified  even  though  the  choice  cannot 

be  made  with  confidence  at  present,  especially  in  view  of  the  excessive 

scattering  and  absorption  in  both  single-crystal-  and  polycrystalline- 

fibers.  A  new  extrusion  technique  using  a  quartz  capillary  with  a  square 

exit,  with  hydrogen-chloride  flushing  and  argon  cooling,  is  proposed. 

Intrinsic  Rayleigh  plus  Brillouin  scattering  has  typical  scattering 

-4 

coefficients  at  1  urn  ranging  from  4  x  10  to  0.4  dB/km  for  the  typical 
materials  considered.  Intrinsic  Raman  scattering  can  be  much  greater, 
by  a  factor  of  40  in  the  examples,  than  the  Brillouin  scattering,  which 
is  sometimes  tacitly  assumed  to  set  the  intrinsic  scattering  limit. 
Scattering  from  vacancies  or  impurity  ions,  either  isolated 
or  in  small  clusters,  is  negligible,  but  scattering  from  extended 
imperfections  can  be  great.  Quenching  or  annealing  temporarily 
changes  clustering  and  scattering,  but  clustering  and  water  or  OH" 
uptake  cause  scattering  in  alkali  halides  to  increase  with  crystal  age. 

It  is  shown  that  multiple  scattering  and  scattering  by  the  strain 
field  around  an  impurity  ion  are  negligible.  Laser-damage  thresholds 
are  tabulated. 
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A.  INTRODUCTION 


The  purpose  of  the  program  is  to  perform  theoretical  analyses  of 
fiber  materials  and  fibers  in  order  to  determine  the  practical-  and 
ultimate-intrinsic-limits  of  low-loss  fibers.  Potential  uses  of  low-loss 
fibers  include  long-distance  communications,  power  transmission,  and  such 
devices  as  various  detectors.  The  most  important  limitations  to  attaining 
useful  low-loss  fibers  are  the  losses— by  absorption  and  particularly 
scattering— and  dispersion,  which  limits  the  bandwidth  in  communi cations 
fibers.  The  emphasis  of  the  program  is  on  crystalline  fibers,  not 
necessarily  because  they  were  expected  to  be  superior  to  glass  fibers  ,  but 
because  analysis  of  crystalline  fibers  was  badly  needed. 

The  character  of  the  study  changed  considerably  since  the  study  was 
first  proposed.  Originally,  absorption  and  Brillouin  scattering,  displayed 
in  V  curves  as  discussed  below,  were  discussed  as  intrinsic  limits,  and  the 
practical  limit  was  considered  to  be  the  purification  of  the  materials  to 
reduce  the  absorption.  Since  scattering  in  crystalline  fibers,  both  poly¬ 
crystalline  and  single-crystal ,  is  a  severe  problem  that  must  be  solved 
before  reducing  the  absorption  will  be  useful,  the  scattering  investigation 
became  a  major  part  of  the  program. 

The  loss,  or  attenuation,  is  characterized  by  the  extinction  coefficient 
8ex  *  8ab  +  8sc*  cons^stin9  of  absorption  Sab  and  scattering  B$c  contributions. 
The  lowest  loss  attained  to  date  is  Bgx  =  0.3  dB/km  for  a  fused-silica  fiber 
at  1.6  um.  Other  low  attained  values  of  current  interest  are  2  dB/km  at 
0.85  um  and  0.5  dB/km  at  1.3  urn.  For  commercially  available  fiber,  1.5  dB/km 
for  1  to  1.6  um  operation  is  considered  standard-product  performance.^ 


Recall  that  1  dB/km  =  2.3  x  10  cm-  ,  which  is  an  extremely  low  loss. 

The  value  of  0.3  dB/km  is  believed  to  be  the  intrinsic  1.6  ym  value 
and  the  lowest  attainable  value  at  any  wavelength  for  fused  silica.  After 
many  years  of  research,  the  intrinsic  limit  of  fused  silica  has  been 
reached.  Thus, other  materials  must  be  considered  in  order  to  obtain 
lower  loss.  Both  glasses  and  crystalline  materials  are  candidates.  Crys¬ 
talline  materials  have  both  potentially  lower  scattering  and  absorption, 
but  severe  technical  problems  must  be  overcome  before  the  potential  of 
crystalline  fibers  can  be  realized. 

The  goal  of  the  Clearday  program,  as  the  low-loss  fiber  program  of 

the  Defense  Advance  Research  Projects  Agency  is  called,  was  originally  to 

-4 

reduce  the  loss  to  such  low  values  as  10  dB/km  or  even  lower.  Since 

such  low  values  are  now  considered  to  be  unrealistic  for  the  near  future, 

the  goal  is  to  reduce  the  loss  as  much  as  possible  below  these  current 

values.  For  example,  reduction  in  the  loss  by  a  factor  of  30  would  give 
_2 

Bex  =  10  dB/km  as  a  typical  goal. 

2 

In  selecting  low-loss  materials,  V  curves  are  usually  used.  The 
V  curves  are  plots  of  the  absorption  coefficient  and  the  scattering 
coefficient  as  functions  of  frequency,  as  illustrated  by  the  tabulation 
of  V  curves  in  Sec.  G.  Even  though  the  V  curves  have  been  considered 
extensively  in  the  literature  and  are  useful  for  glasses,  they  are  perhaps 
more  misleading  than  useful  for  crystalline  materials  at  present.  First, 
the  scattering  half  of  the  V  is  the  intrinsic  Brillouin  scattering.  The 
low  theoretical  value  of  the  Brillouin  scattering  in  crystals  was  indeed 
one  of  the  initial  reasons  for  considering  crystals  for  low-loss  fibers. 
However,  as  shown  in  Secs.  B  through  E,  extrinsic  scattering  is  currently 
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several  orders  of  magnitude  greater  than  the  intrinsic  Brillouin  scattering. 
Reduction  of  the  scattering  in  crystalline  fibers,  even  single  crystals,  is  a 
cardinal. problem  that  must  be  solved  in  order  to  attain  low-attenuation 
fibers. 

2  2 

Second,  the  dispersion  (particularly  d  n/dX  ,  where  n  is  the  real  part 
of  the  index  of  refraction  and  X  is  the  wavelength),  rather  than  loss,  often 
limits  the  use  of  fibers  because  signals  carrying  information  necessarily 
have  nonzero  bandwidth  and  the  different  frequencies  within  the  bandwidth 
travel  at  different  velocities. 

2  2 

The  importance  of  the  material  dispersion,  as  d  n/dx  will  be  called, 
can  be  seen  as  follows:  The  group  velocity  is 

,,  _  dw  _  c  m 

vg  "  dF  “  d6/dkvac  ’ 

where  kvgc  =  2tt/X  ,  X  is  the  wavelength  in  vacuum,  m  =  ckvac,  and  the 

propagation  constant  6  (denoted  by  k  in  some  references)  determines  the 

spatial  phase  (£~exp(i6z)  for  propagation  along  the  z  axis).  Both  the 

propagation  characteristics  of  the  fiber  and  wavelength  dependence  of  the 

refractive  index  n  determine  the  dependence  of  3  on  k  .  The  function 

vac 

B(kvaC)  is  complicated  in  general.  However,  by  considering  the  effects 
of  the  propagation  characteristics  alone  and  the  refractive-index  disper¬ 
sion  alone,  it  can  be  shown  that  the  refractive-index  dispersion  usually 

3 

dominates  the  propagation-charateristic  dispersion. 

The  following  calculation  shows  that  the  refractive-index  dispersion 

2  2 

is  zero  when  the  material  dispersion  d  n/dX  =  0.  In  order  for  the  various 
frequency  components  of  an  informat ion -carrying  laser  signal  to  arrive 
at  the  output-end  of  the  fiber  at  the  same  time,  the  group  velocity  must 
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be  independent  of  the  frequency,  that  is,  dvg/dw  =  0.  For  propagation 
of  a  plane  wave  in  an  infinite  medium  of  dispersive  refractive  index  n(X), 
the  propagation  constant  is 


8  =  nk 


vac 


(2) 


Assuming  an  infinite  medium  is  equivalent  to  formally  neglecting  the 
propagation  characteristics  of  the  fiber.  The  derivative  of  vg  in  Eq.  (1) 
is 


doj 


c  dk 


vac 


(«/dkrac)2 


The  value  of  dvg/du  is  zero  at 

d2(nkvac)/dkvac2  ’  0  <3> 

With  kvac  =  2tt/X.  where  X  is  the  vacuum  wavelength,  it  is  easy  to  show  that 
Eq.  (3)  is  satisfied  when 

d2n/dX2  «  0  (4) 

P  7 

is  satisfied.  The  wavelength  at  which  d  n/dX  =  0  is  called  the  zero- 
dispersion-crossover  wavelength. 

Third,  the  extremely  low  values  of  the  extinction  coefficients  at  the 

minima  of  the  V  curves  for  many  crystalline  materials  surely  cannot  be 

attained  in  practice  in  the  foreseeable  future.  In  addition  to  the  scattering 

limitations  mentioned  above,  the  purity  of  the  crystals  required  to  attain 

the  low  absorption  is  unrealistically  low  in  at  least  most  cases  of  interest. 

-4 

As  discussed  in  Sec.  H,  impurity  concentrations  of  only  10  part  per  billion 
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can  give  rise  to  a  typical  goal  value  of  the  absorption  coefficient, 

_2 

6  =  10  dB/km.  Even  the  measurements  of  both  the  impurity  concentrations 

-4  -2 

of  10  parts  per  billion  and  the  absorption  coefficient  B  =  10  dB/km 

are  beyond  the  current  state  of  the  art. 

In  high-power  laser  windows,  absorption  is  usually  more  important 

than  scattering  because  absorption  leads  directly  to  heating,  which  in  turn 

leads  to  such  failure  as  fracture,  vaporization,  plasma  ignition,  or  optical 

distortion.  By  contrast,  the  total  loss,  including  scattering  as  well  as 

absorption,  is  believed  to  be  the  limiting  factor  in  low-loss  fibers. 

There  are  many  related  topics  that  are  not  considered  in  the  present 

investigation.  These  include:  applications  of  fibers,  refractive-index 

gradation,  claddings,  and  guided-mode  theory. 
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B.  SCATTERING  BY  MACROSCOPIC  IMPERFECTIONS 


In  this  section  it  is  shown  that  scattering  is  the  most  serious 
problem  limiting  the  operation  of  low-loss  fibers  of  crystalline  materials, 
both  polycrystalline  and  single-crystal.  Any  one  of  the  following  extrinsic 
scattering  processes  can  limit  the  total  loss  to  values  orders  of  magnitude 
in  excess  of  the  intrinsic  value: 

•  surface  imperfections, 

•  voi ds , 

•  macroscopic  inclusions, 

•  dislocations, 

•  strains, 

•  grain  boundaries  in  noncubic  crystals, 

•  contaminated  grain  boundaries  or  imperfections,  and 

•  impurities,  especially  in  clusters. 

All  of  these  scattering  mechanisms  are  potentially  more  severe  in  fibers  than 

in  bulk  crystals  because  current  processes  of  forming  the  cyrstalline  fibers, 

such  as  extrusion,  hot  rolling,  direct  growth,  and  others  induce  sources  of 

scattering.  Values  of  the  scattering  coefficients  Bsc  that  are  calculated  in 

the  present  section  are  listed  in  Table  I.  These  very  large  values  of  S  — 

-5 

ranging  from  -0. 01  to  -100  dB/m--compared  with  the  typical  goal  of  10  dB/m, 
clearly  illustrate  the  severity  of  the  scattering  problem. 

A  complete  formalism  of  the  scattering  problem  is  given  in  Sec.  E.  The 
simple  estimates  in  the  present  section  suffice  for  current  use  because  the 
scattering  is  so  strong  that  the  simple  estimates  are  sufficiently  accurate 
to  show  that  the  scattering  is  overpowering. 
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First  consider  scattering  by  voids  or  macroscopic  inclusions.  The 

simple  model  of  spherical  particles  of  refractive  index  np  and  density  N 

particles  per  unit  volume  imbedded  in  a  host  dielectric  of  refractive  index 

is  used.  The  absorption  coefficient  for  this  Rayliegh-Mie  scattering  is 
12  3 

well  known.  '  *  A  typical  scattering  cross  section  os  of  a  single  particle  of 

radius  a  as  a  function  of  kHa,  where  kH  =  nHk  with  k  =  2tt/X,  is  sketched  in 
Fig.  1. 

The  particles  of  interest  have  ka>l;  thus  the  scattering  cross  section 

.  2 
is  approximately  equal  to  twice  the  geometrical  cross  section  .  As 

an  example  of  ka  >  1  being  satisfied,  solving  k^a  =  1  for  a  gives 

a  =  X/2irnH  (1 ) 

with  nH  *  1.5,  Eq.  (1)  gives  a  =  0.11X,  and  for  X  =  1  ym,  the  value  of  a  is 
0.11  pm.  Particles  with  a  <"  X/2-rrn^  have  smaller  cross  sections,  as  seen  in 
Fig.  1  and  discussed  below. 

The  scattering  coefficient  Bsc  is 

esc  =  °SN  <2a) 

■’  2ira2N  ,  for  ka  ">  1  (2b) 

=  (4/27)  (nr^  -  n^)^  (k^a)4  ua^N,  for  ka  <  1  (2c) 

Equation  (2a)  is  the  standard  expression  for  3sc  in  terms  of  a$;  Eq.  (2a) 

~  2  12  3 

follows  by  setting  og  =  2Tra  ;and  Eq.  (2c)  is  the  usual  Rayliegh  result.  ’  ’ 

1  2 

The  factor  of  two  in  the  cross  section  is  discussed  in  standard  texts.  ’ 

For  a  typical  case  of  1  pm  -  radius  particles  or  voids  spaced  200  pm  apart, 
the  density  is 
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N  =  (1/200  pm)3  =  1.2  x  io5  cm'3, 

and  for  ka  >"  1  (that  is,  X  ~  9.4  urn  for  nH  =  1.5),  Eq  (2b)  gives, 

Bsc  =  7.8  x  io-3  cm'1  =3.4  dB/m 
5  -3 

Thus,  the  small  concentration  of  10  cm  of  1  pm  radius  particles  gives 

rise  to  a  scattering  coefficient  that  is  fi«  orders  of  magnitude  greater 

-2  5  -3 

than  a  typical  goal  of  10  dB/km.  The  concentration  1.2  x  10  cm  and 

radius  a  =  1  pm  correspond  to  a  fractional  volume  fy  of 

fV  ■  Vticles'Vsample  * 

=  5.2  x  10^  =  0.52  ppm. 

This  small  fractional  volume  corresponds  to  the  very  large  value  of  3.4  dB/m. 

The  fiber  must  be  essentially  particle  and  void  free  in  order  to  attain  the 
goal  of  B$c  «  1  dB/km. 

In  the  calculation  above,  it  was  tacitly  assumed  that  the  scattering 
cross  section  asc corresponds  to  scattering  out  of  the  fiber.  This  approximation 
is  well  satisfied  for  the  purpose  of  the  estimate  because  the  scattering 
angle  for  a  single  sphere  is  approximately  Trn^,  where  n&  =  jn^  -  nr| ,  as  will 
be  shown  below.  This  approximate  value  gives  a  scattering  angle  of  tt/2  for 
nH  =  1.5  and  nr  =  1.0.  Even  though  the  approximation  breaks  down  for  such 
a  large  scattering  angle  (because  the  scattering  angles  at  the  front  and 
rear  surfaces  of  the  sphere  were  assumed  to  be  equal),  the  estimate  indicates 
that  a  single  scattering  is  sufficient  to  scatter  most  of  the  rays  outside 
the  angle  of  total  internal  reflection. 


Next  consider  the  contribution  to  the  scattering  coefficient  from 

surface  imperfections  on  the  fiber.  Somewhat  greater  care  is  required  to 

obtain  reliable  estimates  for  the  surface-induced  scattering  than  for  the 

scattering  by  voids  or  inclusions  in  the  bulk  of  the  fiber.  Nevertheless, 

in  view  of  the  overwhelming  effect  of  scattering  in  the  spherical -parti cal 

model  above,  it  is  intuitively  obvious  that  the  fiber  surfaces  must  be 

essentially  free  of  imperfections  such  as  the  "fish-scales"  observed  by 

4 

Harrington  and  coworkers  and  voids  or  protrusions  on  the  fiber  surface. 

5 ,( 

Single-mode  fiber  operation  is  of  interest  in  communications.  Marcuse 
has  shown  that  surface  imperfections  are  a  strong  source  scattering  loss  in 
single-mode  operation  of  both  round  dielectric  waveguides,  including  fibers, 
and  dielectric-slab  waveguides.  For  example,  "the  radiation  losses  caused 
by  surface  roughness  of  a  fiber  designed  for  single  mode  operation  at  1  urn 
wavelength  can  be  as  high  as  10  dB/km  for  an  ms  variation  of  the  fiber 

o 

wall  of  as  little  as  8  A."  This  and  other  examples  obtained  for  special  cases 
of  Marcuse's  general  results  (as  in  Table  I)  clearly  indicate  that  fiber 
surfaces  must  be  nearly  perfect  in  order  to  attain  even  the  state-of-the-art 
values  of  0.3  dB/km. 

At  the  risk  of  oversimplification,  the  following  simplified  estimate 
explains  the  exceptionally  strong  scattering  from  large-diameter  multimode 
scaley  fibers  fortuitously  well.  A  typical  scattering  coefficient  for  a 
sea ley  fiber  is 

8$c  =  20  dB/m  (3) 

Q 

for  a  250  ym-radius  multimode  potassium-chloride  fiber. 

The  strongly  multimode-fiber  irradiance  is  considered  as  a  collection 
of  rays  that  can  be  treated  by  ray  tracing,  with  diffraction  added  formally. 
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A  ray  through  the  center  of  a  circular  fiber  making  a  small  angle  0  with 
the  fiber  axis  travels  a  transverse  distance  2a^,  where  af  is  the  fiber 
radius,  in  a  distance  L,  where 

0L  =  2af  .  (4) 

The  number  of  times  that  the  ray  reflects  off  the  surface  before  the 

irradiance  is  reduced  by  a  factor  of  1/e  is  denoted  nD.  For  surface 

b 

imperfections  of  dimension  a  (the  radius  of  a  hemispherical  pit,  for  example) 

sufficiently  large  for  kHa  ^  1  to  be  satisfied,  the  large-scattering  cross 
2 

section  os  =  2ira  for  a  spherical  void  in  an  infinite  medium  suggests  that 

the  scattering  is  strong  at  the  surface.  Thus,  a  reasonable  value  of  nD  is 

5 

nB  =  4,  for  kHa  >  1 .  (5) 

Formally  reducing  the  scattering  by  a  factor  of  (4/27)(nH2  -  l)2(kHa)4 
-  e  for  k^a  <  1,  according  to  Eq.  (2c),  gives  the  following  equation  for 

nB 

nR/4 

(1  -  e)  8  =  1/e  (6) 

For  e  «  1 ,  Eq.  (6)  gives 

nB  *  4/e  =  4[(4/27)(nH2  -  l)2^)4]"1  (7) 

Setting  Bsc  =  (nBL)  ^  and  using  Eqs.  (4),  (5),  and  (7)  gives 

Ssc  =  0/8a^  ,  for  k^a  ">  1  (8a) 

=  (1/54) (nH2  -  l)2(kHa)4  e/af,  for  kHa  <  1  (8b) 
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The  value  of  6  to  be  used  in  Eq.  (8)  is  determined  by  the  divergance 
and  alignment  of  the  incident  laser  beam,  the  deviation  from  straightness 
of  the  fiber,  and  diffraction.  For  the  diffraction  limit,  the  value  of 
8  is  formally 


9  =  1.2X/2af 

(9) 

For  the  case  of  a^  =  250  pm,  X  =  10.6  pm,  and  kHa  >  1,  Eqs. 

give 

(8a)  and  (9) 

9  =  1.27  x  10"2  rad 

(10) 

6SC  =  0.254  cm"1  =  55  dB/m 

(ID 

in  excellent  agreement  with  the  typical  experiment  value  of  8  =  20  dB/m 

in  Eq.  (3). 

Next  consider  the  optical  scattering  by  dislocations.  The  order  of 
magnitude  of  Bsc  is  obtained  by  the  following  simple  estimate  of  the  scattering 
by  the  stress  field  associated  with  a  dislocation.  The  magnitude  of  the 

9 

stress  in  the  vicinity  of  an  edge  dislocation  is 


.  _  _  _  Gb 

°pp  '  aee  '  '  2  TT“7 

Gb  cos9 

aP6  =  2  (T  -")  T~ 


(12) 

(13) 


where  G  is  the  shear  modulus,  b  is  the  Burgers- vector  amplitude  of  the 
displacement,  v  is  the  Poisson  ratio,  and  p  and  0  are  cylindrical  coordinates. 

For  the  order-of-magnitude  estimate,  the  angle  factors  in  Eqs.  (12) 
and  (13)  are  formally  neglected,  giving 


I 


<  i 
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! 


a  =  Gb/27r(1  -  v)p, 

=1.4acb/p  (14) 

where  we  used  the  typical  value  of  v  =  0.3  in  the  second  equality.  The 
stress  reaches  the  critical  value  a.  =  G/2tt  for  which  the  lattice  becomes 
unstable  at  the  radius  p  b.  Thus,  the  model  for  the  stress  is 
basically  valid  at  distances  from  the  dislocation  that  are  many  times  the 
lattice  spacing  b,  which  is  not  surprising  in  view  of  the  severe  lattice 
disruption  at  the  dislocation  itself. 

As  a  further  approximation,  the  stress  field  is  approximated  by  an 
infinitely  long  cylinder  of  uniform  stress.  The  radius  a  of  the  cylinder 
is  chosen  as  the  value  p  at  which  o  =  ac/10.  With  this  value  of  a, 

O 

Eq.  (14)  gives,  with  typical  values  of  y  =  0.3  and  b  =  3  A, 

a  =  10b/(l  -  v)  =  4  x  10'7  cm  . 

The  average  stress  in  the  cylinder  is  approximately 

<  a  >  *  dpp[acb/(l  -  v)o]/f*dpp  =  2boc/(l  -  v)a 
=  ac/5 

The  index  of  refraction 

nr  ■  1  +  q  <a>  *  1  +  qcrc/5  (15) 

then  has  a  typical  value 

nr  ■  1  +  3  x  io"3 


I 
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for  a  typical  value  of  ac  =  G/2it  =  3  *  1011(dyne/cm2)/2Tr  (=  30  GPa/2iT  = 

4.4  x  107  psi /2tt )  and  a  typical  value  of  the  strain-optic  coefficient 
q  =  3  x  10“13  cin2/dyne  (=  3  *  lO'^Pa"1). 


The  scattering  cross  section  per  unit  length  of  a  thin  cylinder  of 

2 

2 


index  nr  and  radius  a,  for  ka  ?  1,  is 


C  =  2a 


TT2(ka)3 


(16) 


The  scattering  coefficient  is' 


Ssc  •  N„C  (17) 

where  is  the  number  of  cylinders  per  unit  area  (the  dislocation  density). 
With  X  =  1 .06  ym ,  nr2  +  1  ^  2,  and 

nr2  -  1  =  (1  +  qoc/5)2  -  1 

=  2qac/S  =  6  x  IQ-3,  (18) 


Eq.  (16)  gives 


C 


2(4  x  10“7)fr2 
4 


■  2.4  x  10-16  cm  . 


(19) 


Since  fibers  are  filled  with  a  large  number  of  dislocations  generated  by 
the  manufacturing  process,  we  choose  the  large  value 
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Nd  =  1012  cm"2  (20) 

Substituting  Eqs.  (19)  and  (20)  into  Eq.  (17)  gives 

Bsc  =  2,4  *  10"4  cm_1  =  103  dB/km 

This  estimated  value  of  B  is  sensitive  to  changes  in  the  model, 

sc 

particularly  the  choice  of  the  values  of  the  parameters.  Typical  dislocation 
densities  Nd  range  from8  102  to  10^2  cm-3.  The  choice  of  the  value  of  "a" 

4 

is  also  rather  sensitive  because  C  ~  a  .  Finally,  at  the  longer  wavelength 
10.6  urn,  with  C  -  X'3,  B$c  =  0.10  dB/km. 

Next  consider  the  scattering  by  the  large  inhomogeneous  strains  induced 
in  crystalline  materials  by  the  fiber  manufacturing  process.  The  strains  could 
be  related  to  grains,  but  are  not  directly  related  to  the  grains  in  general. 

The  inhomogeneities  are  modeled  by  a  system  of  spheres  of  constant  stress. 

The  scattering  of  an  individual  sphere  is  estimated,  then  the  total  scattering 
angle  is  estimated  as  the  random-walk  value  of  N$1/2  9^  where  N$  is  the 
number  of  spheres  and  9-j  is  the  average  scattering  angle  for  one  sphere. 

When  the  total  scattering  angle  is  sufficiently  large  that  the  ray  is  no 
longer  totally  internally  reflected,  the  ray  is  assumed  to  be  lost  by  trans¬ 
mission  out  of  the  fiber. 

Consider  the  scattering  by  a  sphere  of  radius  a  and  constant  stress  o. 

The  angle  of  incidence  0^  for  a  ray  having  impact  parameter  l  is 

9.  =  sin"^£/a.  (21a) 

From  Snell's  law,  the  angle  of  refraction  er  is  given  by 
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0r  =  sin'1  (nr  sin  ei)  =  sin'1  (nrt/a) 


The  deflection  angle  is 


0,  =  04  -  0*. 

dir 


(21b) 


(21c) 


Since  the  ray  undergoes  a  second  deflection  upon  leaving  the  sphere,  the 
total  scattering  angle  0-j  for  6^  <<  1  is 


With 


Qt  =  2(0.  -  0r) 


nr  =  1  +  V 


Eqs .  (21)  give 


(22) 


0-,  =  2 1 sin-1  {1/ a)  -  sin"1  {1/ a  +  nA£/a)| 

^  2|sin-1  (1/ a)  -  sin'1  {1/ a)  -  nAU/a)(l  -  l2/a2)~V2 
=  2nA£(a2  -  £2)"1/2 

For  uniform  illumination  of  the  sphere,  the  probability  of  a  ray  having 
impact  parameter  between  l  and  l  +  is 

2tt££.  a. 

A  9  A 

— r  =  2  x 

Tra  a 


Thus,  the  expected  scattering  angle  is 
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H$  =  7.2  x  1C‘2/nA2  (27) 

For  a  distance  4a  per  scattering,  the  scattering  length  is 

3  "1  =  4aN  . 

sc  s 

With  Eq.  (27),  this  gives 

6SC  ■  3'5nA2/l  •  l28) 

The  change  in  the  index  of  refraction  in  Eq.  (22)  is  obtained  from 
the  approximation 


n  =  1  +  qo 


Comparison  with  Eq.  (22)  gives 

nA  =  qo  .  (29) 

Substituting  Eq.  (29)  into  Eq.  (28)  gives 

Ssc  =  3.5q2o2/a  ■  (30) 

For  the  typical  value  q  =  3  *  10'12Pa_1  and  a  very  large  value  of 

a  =  a  =  5  x  io9Pa  the  value  of  n.  is  1.5  x  10-2,  and  Eq.  (30),  with 
c  a 

a  =  af/2  =  125  pm,  gives 

3  =  3.5(1 .5  x  10_2)2/1 25  x  10’4  cm 

sc 

=  6  x  10"2  cm-1  =  30  dB/m  (31) 

This  large  value  of  8  =30  dB/m  is  attained  for  a  fiber -that  is  severly 

SC  a 

strained— to  the  limit  of  the  theoretical  strength.  For  a  strain  a  =  10  psi 

=  6.9  x  l07Pa,  the  corresponding  value  of  3CS  is 
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Bsc  »  400  dB/kM  (32) 

_2 

As  a  final  example,  a  value  of  Bsc  =  10  dB/kM  corresponds  to  a  = 

3 

0.024  psi  =  1.7  *  10  Pa,  which  is  an  extremely  small  value.  The  examples 
illustrate  that  the  fibers  must  be  essentially  strain  free  in  order  to 
attain  the  goal  values  of  (3$c  «  0.1  dB/km. 

Next  consider  scattering  by  crystalline  grains.  For  the  first  case  of 
noncubic  crystals,  a  ray  at  non-normal  incidence  on  a  grain  boundary  is 
deflected  because  the  index  of  refraction  for  the  ray  in  the  two  grains  is 
different  because  the  direction  of  propagation  with  respect  to  the  crystal 
axes  is  different.  The  scattering  angle  at  a  grain  boundary  is  of  the 
order  of  the  difference  n  in  the  refractive  indices  of  the  two  grains. 

For  example,  Snell's  law  is 

nr  sin  e.  *  (nf  +  n&)  sin  (9.  -  ed) 

=  (nf  +  nA)[sin  0.  -  9d  d  sin  6i.d9i  +  . . .] 

Equating  the  zeroth-order  terms  gives  the  identity  nr  sin  9^  =  np  sin  9.., 
and  equating  the  first-order  terms  and  solving  for  9^  gives 

8d  =  nA  tan  9i  (33a) 

=  n^  ,  at  0.  =  */4  .  (33b) 

As  shown  in  the  similar  calculation  above,  the  average  of  0-j  =  20^  over  a 
sphere  is 

<el>  =  ™a 


(34) 
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Also,  6SC  is  given  in  terms  of  <6^>  by  Eq.  (24),  and  3  is  given  by  Eq.  (28). 

-2 

With  n^  =10  ,  as  a  typical  value  of  the  difference  between  the  indices 

for  propogation  along  different  directions  in  a  noncubic  crystal,  and  a  = 
a^/2  =125  pm,  Eq.  (28)  gives 

S$c  =  2.8  x  10'3  cm"1  =  1.2  dB/m  (35) 

This  large  value  of  8  clearly  demonstrates  that  noncubic  polycrystalline 

J  V* 

fibers  will  have  excessive  scattering.  They  are  not  useful  for  low-loss  fibers. 

Next  consider  the  scattering  by  crystalline  grains  in  cubic  crystals. 
Theoretically,  perfect  grain  boundaries  in  cubic  crystals  do  not  cause  scat¬ 
tering  because  the  value  of  the  index  of  refraction  is  independent  of  the 
angles  between  the  wavevector  of  the  radiation  and  the  crystal  axes.  Con¬ 
taminated  grain  boundaries  can  of  course  cause  scattering,  as  discussed 
in  Secs.  C  and  D.  Also,  a  constant  strain  can  cause  scattering  at  grain 
boundaries  even  in  cubic  crystals  because  cubic  crystals  are  anisotropic 
in  the  strain-optic  coupling.  A  constant  stress  on  a  polycrystalline  sample 
of  cubic  material  can  also  give  rise  to  inhomogeneous  strains  as  a  result 
of  the  gain  structure  and  the  anisotropic  elastic  constants. 

Scattering  by  impurities,  clusters  of  impurities,  voids,  and  contaminated 
grain  boundaries  or  other  contaminated  imperfections  are  considered  in  Secs. 

C  and  D. 


Table  I.  Tabulation  of 


sc 

(dB/m) 


voids  or  inclusions 


1.7 


surface  imperfections 
Marcuse  single  mode 


0.01 


2.7 


surface  imperfections 
multimode 


55 


dislocations 


103 


strains 


30 


grain  boundaries  in 
noncubic  materials 


1.2 


typical  goal 


10 


-5 


coefficients  estimated  in  Sec.  B 


28 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  illustration  of  a  typical  extinction  cross  section 
for  scattering  by  a  sphere. 


I 


v*-.' 


Normalized  scattering  cross  section,  0  /71 ra^  (dimensionless) 
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Figure  1.  Schematic  illustration  of  a  typical  extinction 
cross  section  for  scattering  by  a  sphere. 


Normalized  wavevector,  ka  (dimensionless) 
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C.  DEFECT  CHEMISTRY  OF  ALKALI  HALIDES  AND  THEIR  OPTICAL  PROPERTIES 

W.  J.  Fredericks 

Department  of  Chemistry 
Oregon  State  University 
Corvallis,  Oregon  97331 


Intrinsic-Brillouin  and  intrinsic-Rayleigh  scattering,  which 
determine  the  lowest  possible  scattering  that  can  be  attained  in  fiber 
materials,  have  typical  scattering  coefficients  £sc  =  %  +  ^-jnt  R 
of  1  x  10~9  cm-^  =  4  x  io"4  dB/km  for  lithium  flouride  to  B  x  10~7  cm'^  = 

0.4  dB/km  for  zinc  selenide,  both  at  wavelength  X  =  1  ym.  The  total  intrinsic 
scattering  coefficient,  for  Brillouin  and  intrinsic  Raman  scattering, 
can  be  much  greater--by  a  factor  of  40  in  the  examples  considered--than 
the  Brillouin  scattering,  in  contrast  to  the  tacit  assumption  made  in  the 
the  standard  practice  of  using  the  Brillouin  scattering  in  the 
V  curves  for  crystalline  materials.  Electronic  scattering  from 
vacancies  or  impurity  ions  that  are  either  isolated  or  in  small 
clusters  is  negligible,  but  electronic  scattering  from  extended 
imperfections,  such  as  precipitates,  grain  boundaries,  Suzuki  phases, 
and  dislocations  can  be  great.  The  substantial  variations  in  both 
the  laser  breakdown  thresholds  and  the  amount  of  scattering  that  is 
observed  as  various  regions  of  a  crystal  are  probed  are  explained  in 
terms  of  the  distributions  of  impurities  and  imperfections--which 
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change  in  time.  Quenching  temporarily  reduces  clustering,  thereby 
reducing  scattering  and  possibly  isolated-spot  laser  damage;  whereas 
annealing  increases  clustering,  thereby  increasing  scattering  and 
possibly  isolated-spot  laser  damage.  Clusters  are  expected  to  form 
in  time  even  in  quenched  crystals,  and  0H“  and  substantially 
penetrate  alkali -halide  crystals  if  they  contain  a  few  parts  per 
million  of  divalent  ions  (which  are  used  to  strengthen  alkali  halides). 
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I.  INTRODUCTION 


The  purpose  of  this  report  is  to  consider  the  effect  of  complex 
defects  present  in  alkali  halides  on  their  optical  properties. 

In  spectral  regions  far  removed  from  any  optical -absorption  band, the 
limiting  transparency  of  a  substance  is  determined  by  the  amount  of 
light  it  scatters.  Light  is  scattered  from  any  region  that  varies 
in  polarizability  from  the  neighboring  region  of  the  substance.  Both 
intrinsic  scattering  and  extrinsic  scattering  are  considered.  Sources 
of  extrinsic  scattering  are  discussed.  Significant  results  are  obtained, 
as  summarized  in  Sec.  V. 

The  sections  of  this  report  are  as  follows: 

I.  Introduction 

II.  Intrinsic  Scattering 

III.  Defect  Chemistry  of  Alkali  Halides 

IV.  Defect  Scattering 
Summary  and  Discussion 


V. 
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II.  INTRINSIC  SCATTERING 


In  this  section  it  is  shown:  that  even  intrinsic  scattering 
can  limit  the  loss  of  fibers;  that  the  scattering  in  even  the  best 
current  bulk  materials  is  greater  than  the  intrinsic  value;  and  that 
intrinsic  Rayleigh  scattering  can  be  much  greater  than  intrinsic 
Brillouin  scattering  in  materials  other  than  alkali  halides. 

The  limiting  case  of  low  scattering  occurs  when  only  the  thermal 
fluctuations  in  density  cause  scattering  of  light.  Landau  and  Placzek1 
showed  that  the  intensity  of  the  intri nsic  wavelength-unshifted 
Rayleigh-scattering  light  to  the  intensity  of  the  intrinsic  shifted 
Bril louin-scattered  light  was  a  constant,  ri[_p-  This  quantity,  known  as 
the  Landau-Placzek  ratio,  is  usually  expressed  as 


^LP  Ir/2Ib  =  Y  -  1  =  (By  "  >  (1) 


where  y  is  the  ratio  of  the  specific  heat  capacity  at  constant  pressure 

to  that  at  constant  volume  and  6^  and  6^  are  the  isothermal  and 

2 

adiabatic  compressibilities,  respectively.  Munster  provides  a 

3 

concise  review  of  early  development  and  also  discusses  later  work  . 

Light  scattering  has  been  used  to  study  many  theoretical  and  practical 

problems,  mostly  in  amorphous  materials,  which  have  been  extensively 

•  .4,5,6 

reviewed  ’  . 

Most  discussions  of  scattering  in  crystalline  solids  have  concerned 


7  8 

the  hypothetical  "pure,  perfect  crystal.  *  "  For  the  present  discussion, 

that  case  represents  the  ideal  limit;  all  real  crystals  will  exceed 

0 

this  limit.  Landau  and  Lifshitz  give  a  convenient  expression  for 
nLP  as 


nLP  =  C11  T  (3aL)/cp  {2) 

where  c^  is  a  component  of  the  elastic  tensor,  c*L  is  the  coefficient 
for  thermal  expansion  and  they  assumed  the  index  of  refraction,  n, 
is  independent  of  temperature,  T,  at  constant  density,  p,  (i.e.  (3n/;-T)^ 
=  0).  A  more  accurate  expression  has  been  given  by  Wehner  and  Klein^ 
as  a  modification  of  Landau's  and  Lifshitz1 s  formula  which  considers 
the  non-zero  value  of  (3n/3T)p  and  from  the  fact  that  the  transverse 
strain  must  vanish  in  a  plane  longitudinal  wave.  Instead  a  transverse 
stress  is  present.  In  their  analysis  they  reduce  these  considerations 
to  a  dimensionless  number,  r,  which  can  be  evaluated  from  the  elastic 
constants,  the  linear  thermal  expansion  coefficient,  the  elasto-optic 
coefficients,  (3c/3u)j,  and  the  thermo-optic  coefficients,  (3e/3T)u , 
where  u  is  a  displacement.  They  show  that  the  corrected  scattering 
ratio  nWK  is  given  as 


nWK  =  nLP  0 


(3) 


where  r  can  be  expressed  as 
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when  the  coefficients  and  (3e/9T)p  are  for  hydrostatic  conditions  and 
the  fluctuations  are  uniaxial.  The  two  thermo-optic  coefficients  can 
then  be  related  by 


(3e/3T)u  =  (3c/DT)  + 


+  3Pio) 


where  e. ,  p..  are  components  of  the  dielectric  constant  and  the 

I  I  J 

elastic-optic  (Pockels)  coefficient  tensors,  respectively .  Alternately, 
r  can  be  evaluated  from  the  index  of  refraction,  n,  if  the  appropriate 
coefficient  are  available.  The  corresponding  relations  are 


k i  **c 1 2^cn 


UyaT) 

( 9  n^j  /9  U2  /  y 


and  the  density-independent  thermo-optic  coefficient  can  be  related 
to  the  pressure  independent  coefficient  by 


(7) 


when  coefficients  are  for  same  conditions  as  in  Eq.  (5). 

Wehner  and  Klein  have  evaluated  r  at  room  temperature  for  a  variety 
of  crystals.  In  Table  I  the  ratio  has  been  calculated  for  a  few 
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selected  substances.  Most  of  the  data  for  e,  n,  p^j  and  differential 
coefficients  were  for  visible  wavelength  with  the  exception  of  silicon. 

The  magnitude  of  the  Brillouin  scattering  coefficient,  3g,  can  be 
estimated  from  Eq.  (5) 

,  _  Tr2e4kT  i  2(P-i22  +  P442)  p442  1  (Q] 

B  2A4  |  C11  +  C12  +  2C44  C44  | 

using  the  data  in  Table  II  and  the  relations  (p^  -  p-|2)/2  = 
and  (c-jj  -  =  c^4  for  cubic  crystals.  Then  from  the  basic 

definition  of  the  Landau-Placzek  ratio  the  Rayleigh  scattering 
coefficient  can  be  estimated  from  the  relation 

•!r  =  23bV  •  (9» 

The  sum  3  ^  +3^  provides  a  measure  of  the  limiting  loss  of  light 

in  a  substance  through  non-absorption  processes.  At  room  temperature 

Li F  and  NaF  are  the  most  transparent  of  the  substances  considered 

here.  However  at  other  temperatures  or  in  wavelength  regions  where 

the  tails  of  absorption  bands  cause  n  or  c  to  vary  this  may  not  be 

12 

so.  The  latter  problem  has  been  considered  extensively  .  Temperature 
effects  enter  in  a  more  subtle  way.  Pohl  and  Schwarz1^  discuss  the 
exceptional  case  of  NaF.  The  temperature  dependence  of  and  the 
Rayleigh  scattering  coefficient  is  determined  by  p.j,  n,  or  e  and 
(3n/3T)p  or  (3c/jT)  .  (See  Eqs.  (6)  and  (7)).  They  found  aL  to  vary 
from  30  *  10'6  K_1  at  200  K  to  0.01  <  10'6  K'1  at  10  K  and  (--n/^T) 
counteracts  (3n/rn:)j/(0p/r<T)p.  Sodium  fluoride  is  the  only  alkali  halide  in 
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13 

which  this  has  been  reported  .  This  unusual  property  causes  the  Rayleigh 

scattering  to  vanish  when  the  contributions  (3n/3p)x(3p/3T)  and  (3n/3T) 

i  P 

are  equal  and  opposite. 

Table  II  shows  that  for  the  alkali  halides  the  magnitude  of  the 

Brillouin  scattering  exceeds  that  of  the  Rayleigh  scattering.  This  has  not 
8  14  15 

been  observed.  ’  ’  The  ^ost  probable  cause  of  the  failure  to  obtain 

the  expected  result  is  that  these  theories  apply  only  to  pure,  perfect 
crystals  and  no  such  crystals  can  exist.  All  real  crystals,  even  if 
pure,  must  have  some  form  of  defect  structure  to  satisfy  thermodynamic 
requirements.  If  they  contain  impurities, a  rather  complex  collection 
of  impurity-defect  and  impurity-impurity  compounds  can  exist.  Theimer 
and  his  associates18’17’^8’^9  recognized  that  defects  and  impurities 
could  make  a  significant  contribution  to  the  Rayleigh  scattering,  but  at 
that  time  neither  quality  crystals  nor  a  detailed  understanding  of 
defect  interactions  were  available. 
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III.  DEFECT  CHEMISTRY  OF  ALKALI  HALIDES 


In  discussing  the  effectsof  defect-  and  impurity-species  on  light 
scattering,  it  is  convenient  to  retain  the  conventional  subdivisions 
of  intrinsic  and  extrinsic,  but  we  will  also  need  to  distinguish  between 
symmetric  and  unsymmetric  species.  First  we  need  more  general  relations 
for  defect  species. 

a.  Intrinsic  Defects.  A  pure  crystal  can  be  perfect  only  at 
zero  Kelvin.  Any  pure  real  crystal  must  contain  defects  whose  number 
is  determined  solely  by  the  temperature  with  which  the  crystal  is  in 
equilibrium.  Each  chemical  group  of  crystals  will  contain  basic  types 
of  defects  and  their  derivative  species.  The  common  basic  defect  in 
alkali  halides  consists  of  a  vacant  anion  and  vacant  cation  site, 
that  is  a  Schottky  defect.  Each  of  the  vacant  sites  (vacancies)  bear 
an  effective  charge  of  a  sign  such  that  the  vacant  site  attracts 
ions  of  the  kind  that  normally  occupy  it.  Thus, cation-  and  anion- 
vacancies  are  mutually  attracted.  When  they  occupy  adjacent  sites, 
the  Coulomb  energy  is  minimized  and  that  configuration  my  be  slightly 
more  stable  than  two  isolated  vacancies.  Such  a  configuration  of 
oppositely  charged  vacancies  is  called  a  pair.  Usually  an  adjacent 
pair  is  termed  a  nearest-neighbor  pair,  and  when  once  removed 
a  next-nearest-neighbor  pair,  and  so  forth. 

In  addition  to  Coulomb  attraction,  the  energy  of  the  host-defect 
system  is  lowered  by  the  ions  in  sites  adjacent  to  a  vacancy  or  a  pair 
shifting  from  their  normal  sites  because  their  chemical  bonds  to  the 
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missing  ion  no  longer  exist.  This  "relaxation"  extends  the  defect  into 
the  lattice.  Pairs  tend  to  form  clusters  because  they  are  dipolar 
species.  This  lowers  both  the  electrostatic-  and  relaxation-energy 
of  the  system. 

The  various  species  of  intrinsic  defects  can  be  represented  by 
psuedo-chemical  reactions  denoting  the  formation  of  basic  and  compound 
defects.  All  intrinsic  defect  species  can  be  represented  by  the  psuedo- 
chemical  reactions 


Perfect  Crystal  =  V+  +  V~ 
a  c 


(10) 


and 


+  =  [(Va»c)5U  ,  (11) 

where  Va  and  Vc  represent  anion  and  cation  vacancies  respectively. 

Brackets  are  used  to  indicate  a  complex  defect  species.  Here  <5V's 

are  stoichiometry  numbers,  whi le  subscript  6  denotes  the  number  of  vacancy 

pairs  in  a  cluster,  and  t;  indicates  the  type  of  pairs  or  clusters  formed, 

that  is  nearest  neighbor,  next  nearest  neighbor,  etc. 

It  has  been  shown  that  the  concentrations  of  the  various  defects 

20  21 

are  related  through  a  set  of  simple  mass-action  expressions,  ’  which  are 
written  to  include  all  complex  species  as 

a&ac  =  =  exp  Us/k)  exp  (-Az$/kT) 


(12) 
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and 


a$c/aaac  =  KjS  -  zfic  exp  U6c/k)  exp  (-hfic/kT)  .  (13) 

In  these  equations  the  thermodynamic  activity,  enthalpy  and  entropy  of 
formation  are  represented  as  cl,  h  and  a,  respectively.  The  mass-action 
constants  are  given  as  K,  the  subscripts  denote  species  as  given  above 
with  the  additional  subscript  s  to  identify  Schottky- defect  equilibrium. 
The  configuration  entropy  is  z  and  the  Boltzmann  constant  is  k. 

The  experimental  methods  usually  used  to  evaluate  the  thermodynamic 
quantities  do  not  distinguish  between  the  various  compound  defects,  but 
simply  separate  charged  defects  from  neutral  defects  on  the  assumption 
that  all  the  neutral  defects  are  [(V  V  hi],  that  is  nearest-neighbor 

a  C  I 

vacancy  pairs.22  If  species  [(VaVc)-j2j,  [(VaVc)2l]  or  [(VaVc)3l]  are 
stable  in  alkali  halides,  they  have  not  been  detected.  Whatever  species 
are  present  their  activities  are  set  by  the  temperature  at  which  the 
defect  system  was  equilibrated  and  are  given  by 

=  K6K6s  .  (14) 

b.  Extrinsic  Defects.  Even  the  best  purification  methods  leave 
some  foreign  atoms  in  any  real  crystal.  The  defect  species  that  these 
foreign  atoms  form  depends  on  whether  they  incorporate  substitutional ly 
or  interstitially  in  the  host  lattice,  whether  they  are  homovalent  or 
aliovalent.and  their  ionicity  with  respect  to  the  host  lattice. 
Experimental  evidence  indicate  that  all  comnon  impurities  in  alkali 
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22 

halides  are  substitutional ly  incorporated  in  the  crystal.  This 
situation  gives  rise  to  a  relatively  simple  defect  chemistry, which  can 
be  described  by 

+  -  tf>yvc)v)5d  .  (15) 

«?  +  -  [<A^(*,)V)SC]  .  06) 

+  =  [f(M.)E(Ay)v)6]  ,  07) 

+  K  +  ■  [UM#AVKfl<^T"ft}5v-<Vc)5v>iS:i*  .08) 

and 

6M^  +  EV  =  [(M^)6EP]+  .  (19) 

Cations  and  anions  are  represented  by  M  and  A,  respectively,  l  and  j 
denote  impurity  species  while  h  denotes  the  host  species.  With  the 
exception  of  s,  which  specifies  type  of  impurity-complex,  the  lower 
case  Greek  letters  are  stoichiometry  or  charge  numbers.  The  impurity- 
or  defect-charge  is  given  as  the  excess  or  deficiency  of  charge  at 
the  lattice  position.  In  these  equations  the  cluster  size,  6,  has  been 
factored  from  the  stoichiometry  numbers.  When  the  product  of  the 
reaction  forms  a  separate  phase  this  is  indicated  by  the  symbol  4, 
and  when  an  extended  defect,  such  as  a  grain  boundary  or  dislocation,  is 
a  reactant  it  is  indicated  by  EV.  Mass-action  equations  similar  to 
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(12)  and  (13)  can  be  written  for  reactions  (15)  through  (19).  The  most 
significant  characteristic  of  the  defect  chemistry  of  ionic  compounds 
is  that  the  various  intrinsic-and  extrinsic-defect  reactions  are  coupled 
through  sharing  common  defect  of  impurity  species.  This  feature  makes 
it  crucial  that  the  quality  of  the  crystals  used  in  the  study  of  physical- 
or  chemical-properties  dependent  on  defect  species  be  very  carefully 
controlled,  and  contributes  greatly  to  the  remarkable  diversity  in 
the  magnitude  of  thermodynamic  parameters  separated  as  experimental 
evidence  for  impurity- related  properties  in  alkali  halides. 

c.  Characteristic  Impurity  Behavior.  In  alkali  halides, v  seldom 
exceeds  1  because  of  the  energy  required  to  form  the  additional 
vacancies  necessary  to  compensate  the  impurity  charge.  It  has  been 

shown  that  Ce3+  (i.e.  v  =  2)  is  neither  very  soluble  nor  mobile  in 

23  24  3+ 

KC1 .  However,  others  reported  that  Bi  diffuses  with  a 

remarkably  low  migration  energy  and  suggest  an  interstitial  diffusion 

mechanism.  However,  the  DQ  (which  is  proportional  to  the  square  of 

the  jump  distance)  is  inconsistent  with  interstitial  diffusion.  In 

general,  impurities  with  v  >  1  should  not  be  found  in  alkali  halides 

unless  intentionally  added. 

Experimental  measurements  of  t  are  not  sufficiently  sensitive  to 
distiguish  various  values  of  £.  A  common  method  of  assigning  binding 
energies  to  nearest  neighbors  (c  =  0)  and  next-nearest  neighbors  U  =  1) 
for  impurity/Vacancy  complexes  [MVc(t)l  is  from  dipole-relaxation  times 
obtained  from  dielectric  relaxation/0  d.c.  polarization,  and  ionic  thermal 
currents.27  This  work  shows  that  when  the  impurity  radius,  r^  is  larger 


than  or  near  the  radius  of  the  ion  for  which  it  substitutes,  r  ,  the 

w 

predominate  value  of  v  is  0,  but  when  r^  is  much  smaller  than  rn 

significant  concentrations  of  the  complex  [M\>c(l)]  are  formed. 

28 

Recently  Catlow  and  coworkers  have  calculated  the  binding  energies 

for  alkaline  earth-vacancy  complexes  in  KCl ,  NaCl ,  KBr.  They  apply  a 

29 

modern  version  of  the  Mott-Littleton  method  using  a  commercial  program, 
Harwell  Automatic  Defect  Evaluation  System,  commonly  denoted  HADES,  to 
achieve  rapid  convergence  of  the  large  lattice  summations  required.  Their 
results  suggest  that  next-nearest-neighbor  complexes  should  be  much  more 
common  with  large  impurities  than  experimental  work  suggests.  Unfortunately 
they  can  only  calculate  energies  assuming  e  =  h  and  can  make  no  estimate 
of  entropy  so  that  careful  comparison  of  free  energies,  g,  cannot  be 
done  as  the  derivation  of  h  from  the  experimentally  measured  g. 

Of  more  importance  here  is  the  tendency  for  these  complexes 
to  aggregate  into  clusters  as  described  by  Eq.  (16).  Clustering  was 
discovered  in  the  60's.^’^  *32,33,34,35  an(j  QryCjen30  noteC| 

in  alkali  halides  quenched  from  350  to  400°C  all  the  impurity  ions 

are  initially  present  as  dipoles  (two  to  five  percent  dissociated 

into  isolated  M+  and  V~  ions),  then  aggregate  into  trimers  as  the 

crystal  ages.  Trimers  were  reported  to  form  by  third-order  kinetics. 

As  third-order  reactions  are  extremely  rare,  considerable  controversy 
37  38 

has  resulted.  ’  These  kinetic  arguments  are  not  of  immediate 

importance  here,  except  to  note  that  these  aggregation  processes 

can  occur  at  storage  temperatures  for  crystals.  Cappelletti  and 
39 

DeBenedetti  have  found  direct  evidence  from  ionic  thermal  current 
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35 

measurements  that  some  second-order  reaction  occurs.  Crawford  has 
developed  a  model  utilizing  binary  reactions  to  form  trimers  that 
would  appear  to  be  a  third-order  reaction.  The  rate  constant  is 
given  by 

-Sr/k  kT  -(Eu,  -  Er.)kT 

k3(apparent)  =  4ttNl  e  5  Dq  rQ  e  1  e  B  .  (20) 

The  activation  energy  for  the  reaction  as  measured  from  the  temperature 
dependence  of  k3  is  £*  and 

E*=Em-Eb  ,  (El) 

where  Eg  is  binding  energy  for  dimer  formation,  rD  is  the  capture  radius 
of  dimers,  DQ  is  the  temperature  independent  preexpoential  and  EM  the 
migration  energy  for  diffusion  of  monomers. 

Assuming  no  lattice  relaxation  about  the  defect  and  assuming 
that  E  .  the  formation  energy  of  the  impurity-vacancy  complex  is 

a 

1  eV  we  can  easily  estimate  the  relative  electrostatic  energy  of  the 
various  configurations.3^ 

The  trimer  is  postulated30, 3^,34,41 ,51  tQ  ge  six_sidec|  comple 

on  the  (111)  planes  of  the  cation  sub-lattice  with  an  isolated  cation 

35 

in  its  center.  The  total  electrostatic  energy  of  such  an  array  is 
given  by 

2 

43)  =  -  [6  +  |  -  -  ]  »  -4.036  q2/Kr  (21a) 

'  Kr0  /3 
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and 

e|3)  =  e|3)  -  3Ea  =  1.063  e2/<r0 

where  -q  is  the  electronic  charge,  k  the  static  dielectric  constant, 
and  rQ  is  the  cation-cation  distance  in  the  {111}  plane.  For  KC1 
rg  =  4.44  x  ID"8  cm  and  k  =  4.64  which  yields  e|3^  =  0.73  ev  and  for 
NaCl  (rQ  =  3„94  x  10‘8  cm,  k  =  5.62)  e[3^  =  0.668  eV. 

The  dimer, as  a  planar  quadripole  lying  in  the  (100)  planes, has  a 
total  eletrostatic  energy  given  by 

2 

ET2)  =  [2  -  ~  ]  =  2.586  q2/Kr  ,  (22) 

1  Kr0  u 

and 


e[2)  =  e|2)  -  2Ea  =  0.586  c2/KrQ 

The  values  for  KC1  and  NaCl  of  e|2^  are  0.413  and  0.377 eV,  respectively. 

If  the  dimers  form  by  loosely  coupling  as  a  linear  arrangement 
along  <110>  directions  the  pair  of  impurity-vacancy  complexes  may 
be  adjacent  (c;  =  0)  or  next-nearest  neighbor  (r,  *  1)  or  r,  =  3  etc.  Such 
a  pair  of  complexes  has  a  total  electrostatic  energy  given  by 


+  _J _ ? _ 

n  +  2  n  +  1 


+  2], 


(23) 


where  n  is  the  number  of  unit  rg's  the  pair  is  separated  along  the  < 11 0> 
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row.  The  interaction  energy  is,  then,  given  by 


E  =  r!_  rl  +  — ] _ 

tI  <r0  Ln  n  +  2  n  +  1J 


For  n  =  1,  Ej  =  0.235  eV;  n  =  2,  Ej  =  0.059  eV;  n  =  3,  Er  =  0,024  eV 
in  KC1  and  in  NaCI  similar  spacings  give  0.214  eV,  0.054  eV  and  0.021  eV, 
respectively. 

On  the  basis  of  this  crude  estimate  if  we  neglect  the  contribution 
of  the  entropy  to  the  free  energy,  the  concentration  of  hexipoles  and 
quadripoles  should  not  differ  greatly  in  concentration  in  the  temperature 
range  of  50-1 00°C. 

Manganese  in  sodium  chloride  is  the  most  extensively  studied 
system,  followed  by  cadmium  in  NaCI.  ’  ’ 

Unfortunately  there  is  little  agreement  on  the  enthalpies  or  entropies 
derived  from  the  various  experiments,  or  even  universal  agreement  on  the 
species  that  exist.  Most  workers  accept  trimers  as  the  stable  species 
and  assume  very  rapid  formation  of  dimers  followed  by  a  slightly 
slower  buildup  of  trimers.  However,  Lil  ley^’59’^^4  puts  forth 
strong  arguments  that  cast  doubt  on  their  existence.  Our  concern 
is  only  the  fraction  of  typical  impurities  found  in  aggregates  or 
in  large  defects  and  their  rate  of  formation.  One  wouTd  expect  that 
with  the  large  amount  of  work  done  on  such  systems,  concentrations  could 
be  calculated  with  some  confidence.  Unfortunately  that  is  not  true. 

In  Table  IV  a  selection  of  thermodynamic  parameters  are  presented  for 
a  few  impurity  systems.  These  large  species  must  be  formed  by 
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diffusion  of  the  divalent  ion  and  the  vacancy  to  adjacent  sites. 

This  occurs  with  remarkable  speed.  Table  V  gives  frequency  factors 
and  activation  energies  for  formation  rates  of  the  species  assumed 
to  be  {MVc(0))3. 

The  salient  feature  of  these  large  impurity  species  is  the 

oo  cc 

remarkable  speed  with  which  they  form  at  low  temperatures.  Symmons  ’ 

-4  ?+ 

notes  that  in  NaCl  containing  less  than  10  molar  fraction  of  Mn  (below 

2+ 

the  concentration  at  which  Mn  forms  visible  precipitates),  the 

equilibrium  between  (MnVc(O)},  and  (MnVc(l)},  was  established  in  15  days 

at  -120°C.  At  40°C  equilibrium  would  be  achieved  in  70  seconds. 

2+ 

Symmons  finds  less  than  0.02  of  the  Mn  in  NaCl  present  as  single 

dipoles  when  the  system  has  reached  equilibrium  at  21 °C.  Cappelletti 
27 

and  Fieschi  claim  that  it  requires  100  hours  for  equilibrium  with 
"microprecipitates"  in  KBr:  Sr  at  61°C.  To  achieve  equilibrium  between 
nearest  neighbor  and  next-nearest  neighbor  species  requiresa  host  cation  to 
exchange  positions  by  one  jump.  The  jump  frequencies  from  nearest  neighbor 
and  next-nearest  neighbor  is  given  by 


4  =  1.059  x  1014  e-°-6711/kT 


(24) 


and  the  jump  from  next-nearest  neighbor  to  nearest  neighbor  is  given  as 


W3  =  7.357  x  1013  e-°-6319/kT 

Clearly,  a  one-jump  process  can  achieve  equilibrium  rapidly. 

Species  such  as  trimers  or  microprecipitates  that  require  the 
migration  of  impurity  ions  into  a  cluster  form  in  much  shorter 
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times  than  would  be  expected  if  the  impurities  are  reasonably  uniform 

in  their  distribution  in  the  host  and  they  migrate  by  normal  diffusion 
22 

processes.  Unless  a  very  special  effort  is  made.it  is  almost  impossible 

to  grow  an  alkali  halide  with  a  uniform  distribution  of  impurities.  The 

nonuniform  distributions  have  caused  problems  in  many  experiments.**^ 

General  reaction  17  represents  association  reactions  among  cations 

and  anions  within  a  lattice.  These  have  been  observed  to  affect  most 

measurements  made  on  ionic  crystals.  Here  are  a  few  examples  to 

illustrate  this:  ionic  conductivity,70  ionic  thermal  currents,^7 

71  26 

perturbed  optical  absorption,  dielectric  relaxation,  increased 

72  72 

solubility  of  impurities,  and  anomalous  diffusion.  In  addition 

to  direct  reactions,  cation-cation  or  anion-anion  interactions  that 

43 

markedly  affect  the  properties  being  measured  can  occur. 

A  related  process  that  does  not  occur  entirely  within  the  host 

crystal  is  the  modification  of  an  impurity  within  a  crystal  by  a 

73  74  75 

reactive  gas  in  the  external  phase.  ’  ’  Typical  of  this  type 

reaction  is  the  apparent  anomalous  diffusion  in  KC1  of  OH"  in  the 
74 

presence  of  C^.  While  such  reactions  are  known,  no  thermodynamic 
data  is  available  and  only  for  the  system  KC1:  OH:  Hg  have  all  the 
products  been  identified.7^ 

Precipitation  of  an  impurity  to  form  a  metastable  phase  is  represented 
by  reaction  ( 18) .  As  written  £  ion  pairsof  impurity  compound 
precipitate  with  <t>  ion  pairs  of  host  M^A  and  £v  cationic  vacancies 
to  form  a  separate  phase.  The  best  known  examples  are  the  Suzuki 
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phases  which  form  with  some  divalent  impurities  in  the  alkali  halides 
with  the  stoichrometry  [M^A2  *  '  as  a  cu^ic  phase. 

Frequently  these  are  written  as  6  NaCl  •  CdCl 2 »  but  here  it  is  written  as 
a  lattice  anionic  impurity  because  its  structure  requires  two  cationic 
vacancies  in  the  normal  alkali-halide  lattice  which  are  compensated  by 
two  cations  in  the  host  lattice.  This  phase  is  stable  below  500°C  in 
NaCl.78  Many  similar  phases  have  been  identified  in  sodium  chloride: 

NaCl:  SrCl2,  NaCl:  MnCl2>  NaCl:  BaCl2  NaCl:  CaCl2,  LiF:  MgF2.77’78‘79’80’81 

Most  early  work  was  done  by  x-ray  diffraction;  but  now  such  phases  have  been 

82  83  84  85 

observed  by  optical  microscropy,  *  by  gold  decoration  replication ’ 

nr 

and  by  direct  transmission  electron  microscopy.  Most  studies  of 

Suzuki  phases  involve  crystals  containing  the  impurity  ion  in  a 

concentration  far  in  excess  of  its  low -temperature  solubility.  However, 

it  seems  to  be  the  species  that  forms  as  a  precursor  to  precipitation 

of  regions  of  MA2.  The  distributions  of  precipitate  sizes  has  been 

83 

determined  by  replication  measurements  and  is  generally  in  accord 

with  that  predicted  for  a  surface-reaction  limited  Ostwald-ripening 
87 

mechanism.  In  this  mechanism,  the  6NaCl  •  MnCl2  was  best  fit  by  an  F-|(u) 
except  for  the  small  number  of  particles  in  which  R/Rc  was  less  than 
0.5,  where  R  is  the  radius  of  the  particle  considered  and  R^.  is  the 
radius  of  a  particle  which  is  instantaneously  in  equilibrium  with  the 
average  solute  concentration,  that  is  the  radius  of  a  particle  that  is 
neither  growing  nor  shrinking. 

A  carefully  grown  and  well  treated  alkali  halide  can  be  relatively 
free  of  extended  defects  [reaction(19)].  If  such  crystals  are  strained, 
as  in  machining  or  polishing,  dislocations  are  introduced.  Rapid  or 
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rough  growth  from  mechanical  vibrations  or  from  thermal  stresses  will 

also  introduce  significant  numbers  of  dislocations.  When  a  crystal  is 

periodically  deformed, dislocations  sweep  through  it  collecting  vacancies 

88  89  90 

as  they  traverse  it,  thus,  becoming  charged.  ’  ’  The  dislocations 

will  not  only  collect  free  Vc  but  will  remove  associated  vacancies,  i.e. 

decompose  { MVC ( 0 ) } . .  Charges  as  high  as  one  ionic  charge  per  two 

88 

cation  sites  on  the  dislocation  core  have  been  reported.  If  no  jogs 
are  formed,  the  dislocation  can  come  to  equilibrium  only  by  migration 
of  V  at  temperatures  too  low  for  anion  diffusion. 

The  equilibrium  concentration  f  of  cation  vacancies  on  the  core 
will  be  given  approximately  by 

r^T  =  ote^kT  ^9+d  "  H>  ’ 

where  a  is  the  concentration  of  free  cation  vacancies  in  distant  parts 

of  the  crystal,  g+d  is  the  free  energy  of  binding  of  a  cation  vacancy 

to  a  core  site  and  9E/3M  is  the  increase  in  the  electrostatic  energy 

of  the  whole  system  resulting  from  the  addition  of  one  vacancy  to  the 
92  93 

core.  ’  The  charge  cloud  satisfies  the  Debye-Huckel  equation  and 

in  NaCl  has  a  radius  of  a  few  hundred  lattice  spacings.  When  divalent 

impurity  cations  are  mobile  they  may  diffuse  to  the  dislocation  core, 

but  they  can't  move  with  the  dislocation  as  it  glides.  They  act  as 

pinning  points  and  their  accumulation  can  be  studied  by  internal 
94 

friction,  but  not  by  experiments  that  detect  moving  charged 
dislocations.  When  processes  occur  that  produce  free  electrons, 
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such  as  photoionization  of  a  monovalent  impurity,  they  may  collect 

on  dislocations.  Ohkura  found  a  trap  depth  of  0.26  eV  for  a 

91 

dislocation  and  0.48  eV  for  a  jog  in  KBr. 

Another  type  of  extended  defect  in  alkali  halides  is  the  grain 
boundary.  These  macroscopic  defects  are  present  in  poorly  grown 
crystals,  especially  Stockbarger-grown  crystals  that  have  seeded 
poorly.  They  are  introduced  by  gross  deformation  of  crystals,  and 
are  present  in  the  pressed  windows  used  for  preparation  of  infrared  win¬ 
dows.  Grain  boundaries  are  very  active  defects  chemically.  They 

collect  impurities  and  provide  charge  anomalies  in  otherwise  reasonably 

95 

good  crystals.  Whipple  has  developed  the  theory  of  diffusion 
along  grain  boundaries.  Harris96,97 .98,99,100,101  has  ex^ens-jve-)y 

studied  the  properties  of  grain  boundaries  in  alkali  halides.  Diffusion 

95 

measurements  are  in  reasonable  agreement  with  Whipple's  theory.  Diffusion 

can  contribute  significantly  to  the  deterioration  of  the  optical  properties 

of  crystals  containing  grain  boundaries.  Impurity  concentrations  along 

grain  boundaries  far  exceed  the  bulk  impurity  content.^01  At  present,  little 

has  been  reported  on  the  charge  anomaly  that  accompanies  the  grain  boundary, 

98 

but  in  gold- decoration  experiments  Harris  and  Cutmore  found  that 
the  influence  of  the  grain  boundary  could  extend  250  to  300  ym  from 
the  grain  boundary  plane.  Even  grain  boundaries  formed  at  such  small 
angles  that  they  are  difficult  to  detect  produce  a  major  disruption 
in  the  properties  of  the  crystal. 

While  the  major  chemical  properties  are  determined  by  the  long- 
range  electrical  forces, these  alone  do  not  explain  the  behavior  of 
defects  and  impurities  in  alkali  halides  or  other  ionic  crystals. 
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Certainly  one  more  localized  force  is  the  displacement  of  the  host  ions 
about  such  defects.  Typical  of  this  type  of  problem  is  the  volume 

103 

of  formation  of  Schottky  defects  in  alkali  halides.  Faux  and  Lidiard 

calculated  these  volumes  using  an  extended  Mott-Littleton  approach 

and  predicted  that  the  volumes  of  formation  were  less  than  the 

molecular  volumes.  For  example,  in  NaCl ,  -0.69  Vm;  KC1 ,  -0.73  Vm; 

KBr,  -0.51  Vm;  i.e.  the  volumes  of  relaxations  are  negative.  Experimentally, 

both  contractions  and  expansions  have  been  reported.  Early  work 

based  on  density  changes  in  KC1  to  which  CaC^  and  SrCl^  had  been 

104 

added  indicated  DV  was  negative.  However,  recent  experiments 

on  the  pressure  dependence  of  ionic  conductivity  in  both  the  extrinsic 

and  intrinsic  regions10^100,10^’100  gave  volumes  of  formation  greater 

than  Vm,  i.e.  an  expansion  of  the  lattice  around  the  vacancies  due  to 

relaxation  of  the  lattice.  For  simple  substitutional  monovalent 

impurities,  even  those  that  produce  "off-center"  defects  such  as 

Li+  in  KC1  can  be  calculated  with  reasonable  success,109  but  later 

results  have  differed  somewhat  from  experimental  results  for  divalent 

28 

highly  ionic  impurities.  These  seem  to  overestimate  the  fraction 
of  next-nearest  neighbor  complexes  formed  with  impurities  near  the 
size  of  the  host  cation.  Clearly  the  experimental  data  should  be 
reevaluated  using  a  slightly  different  model.  As  yet  no  modern  work  has 
been  done  on  impurities  that  form  only  partially  ionic  bonds  with  the 
host  anion.  Nor  is  there  any  theoretical  explanation  of  the  enhanced 
solubility  of  divalent  cations  in  alkali  halides  containing  0H“  or  other 
anionic  species,  or  of  the  unusual  stability  of  these  impurity  compounds. 


J 


Sec.  C 


54 


IV.  DEFECT  SCATTERING 


Only  Theimer  and  his  associates  have  attempted  to  study  defect 
scattering  from  alkali  halides. 16,17,18,19  ^  t^e  t1-me  0f  ^^eir 
work,  neither  good  crystals  nor  a  very  satisfactory  understanding 
of  the  nature  of  impurities  in  alkali  halides  was  available. 

Apparently  many  of  their  crystals  contained  rather  large  quantities 
of  dislocations,  grain  boundaries,  and  impurities  in  undefined  states 
of  aggregation. 

In  a  crystal  with  a  low  dislocation  content  and  no  grain  boundaries, 
but  with  a  small  content  of  divalent  impurity,  the  defect  species  will 
be  small  compared  to  the  wavelength  of  the  incident  light.  The  defect 
species  act  as  independent  scattering  centers.  At  lew  impurity 
concentrations  the  intensity  of  the  Brillouin  scattering  will  not  be 
greatly  changed  and  provides  a  convenient  reference  so  that  a  modified 
Landau-Placek  ratio  provides  a  useful  experimental  measure  of  defect 
scattering. 

This  scattering  ratio  is  given  by 


8tt2N 

n  “  nWK  „  .4 

BB* 


7  XK  §K  \ 


where  X  is  the  wavelength  of  the  light  incident  on  the  scattering  center 
K,  N  is  the  number  of  ion  pairs  per  unit  volume,  is  the  mole  fraction 
of  defects  of  species  K,  a  is  the  effective  polarizability  tensor  of 

In 

defect  K  in  a  particular  host  lattice  and  §  selects  the  elements  of  5 

r\  K 
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that  contribute  to  the  light  scattered  at  an  angle  0  from  the  incident 
beam  of  specific  polarization  from  each  of  the  allowed  orientations  of 
defect  K  that  are  optically  distinguishable.  Equation  (27)  applies 
to  constant  temperature  measurements  because  xK  depends  on  the  temperature 
at  which  the  defect  system  was  equilibrated  and  the  time  since  the 
system  was  equilibrated.  aK  has  the  values  corresponding  to  the 
temperature  of  the  scattering  measurement, as  does  Bg  and  N.  At  a 
particular  analytical  concentration  C.  of  impurity  i_  the  magnitude 
of  n  will  depend  on  how  the  impurity  is  distributed  among  the  various 
defect  species  K  because  each  has  a  different  value  of  a^. 

Equation  (27)  must  also  contain  a  distribution  function,  p(x,K), 
to  account  for  the  inhomogeneous  distribution  of  defect  species  K 
within  the  crystal.  In  general  p(x.,K)  is  unknown  and  unique  to  each 
crystalline  sample  as  it  depends  on  the  host  substance,  ^ ,  and  the 
thermal  treatment  to  which  the  crystal  has  been  subjected.  Except 
at  phase  boundaries,  such  as  precipitates,  p(x,,K)  must  be  continuous 
with  the  gradients  established  by  the  balance  of  diffusional  forces 
and  the  lowering  of  the  free  energy  of  the  host-impurity  system  by 
forming  aggregate  species.  Scattering  experiments  should  be  useful 
in  studying  the  distribution  of  impurities  in  alkali  halides. 

Many  estimates  of  the  polarizability  aK  assume  it  is  proportional 

to  the  volume  of  the  simple  defect  or  to  the  volume  of  the  ion  removed 

1 9 

from  the  site  which  forms  the  defect.  This  probably  underestimates 
because  it  neglects  the  relaxation  of  the  ions  surrounding  the 
defect.  Of  concern  in  estimating  aK  for  a  particular  defect  is  the 
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extent  of  the  perturbation  of  the  electronic  distribution  around  the 
defect.  This  problem  has  not  yet  been  treated, although  much  effort 
has  been  spent  on  electronic  centers.  If  it  is  assumed  that  the  anomaly  in 
the  polarizability  that  causes  scattering  is  proportional  to  the 
relaxation  volume  around  the  defect,  a  more  realistic  estimate  should 
be  obtained. 

For  sodium  chloride  using  1  nm  light, Eq.  (27)  reduces  to 


n  =  0.244  +  (6  x  1038  cm'6)  xK  9K  cxj^ 

K 


(28) 


For  a  symmetric  defect  such  as  an  isolated  vacancy  or  impurity  ion.g^ 

is  unity.  If  a  20  percent  change  in  n  can  be  measured,  then  for 

symmetric  defects,^  must  be  8  x  10~41  cm6.  Small  defects  such 

as  isolated  vacancies  would  not  contribute  much  scattering  at  normal 

concentrations,  but  the  larger  aggregates  can.  Consider  the  example 

of  one  part  per  million  isolated  vacancies  in  sodium  chloride  or 

another  chloride.  Assuming  that  of  the  vacancy  is  equal  to  that 

2  -24  3  2 

of  the  missing  chlorine  ion  gives  =  (3.7  x  10  cm  )  .  Thus, 

9|<  =  10  6  (3.7  x  10-V  =  1.4  x  10~6^  cm6,  which  is  a  factor 

of  6  x  io^  smaller  than  the  value  of  8  x  10  cm6  required  to  give 


a  20  percent  change  in  n. 
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V.  SUMMARY  AND  DISCUSSION 


Intrinsic  Brillouin-  and  intrinsic  Rayleigh-scattering,  which 
determine  the  lowest  possible  scattering  that  can  be  attained  in 
fiber  materials,  have  typical  scattering  coefficients  Bsc  =  R 

of  1  x  10  ^  cm~^  =  4  x  10“4  dB/km  for  lithium  to  8  x  10-^  cm"^  = 

0.4  dB/km  for  zinc  selenide,  both  at  a  wavelength  X  =  1  ym.  The  ratio 
2r\  of  intrinsic  Raman  scattering  to  the  intrinsic  Brillouin  scattering 
is  non-negligible  in  general,  in  contrast  to  the  tacit  assumption 
made  in  the  standard  practice  of  using  the  Brillouin  scattering 
as  that  of  the  V  curves  for  crystalline  materials.  The  total  intrinsic 
scattering  coefficient  Bg  +  Bint  R  can  be  much  greater,  by  a  factor 
of  40  in  the  examples  considered, than  the  Brillouin  scattering.  The 
scattering  curves  in  the  V  curves  in  Sec.  G,  which  include  only  Brillouin 
scattering  for  crystalline  materials,  should  be  recalculated  with 
intrinsic  Raman  scattering  included.  The  extrinsic  scattering  in 
fibers  can,  of  course. be  orders  of  magnitude  greater  than  the  intrinsic 
scattering  as  discussed  in  Sec.  B. 

Electronic  scattering  from  vacancies  or  impurity  ions  that  are 
either  isolated  or  in  small  clusters  is  negligible,  but  electronic 
scattering  from  extended  imperfections,  such  as  precipitates,  grain 
boundaries,  Suzuki  phases,  and  dislocations  can  be  great. 

The  substantial  variations  in  both  the  laser  breakdown  thresholds 
and  the  amount  of  scattering  that  is  observed  as  various  regions  of 
crystal  are  probed  can  be  understood  as  follows:  Even  though  the 
brief  review  in  Sec.  Ill  of  the  current  state  of  impurity  chemistry 
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in  the  simplest  ionic  solids,  the  alkali  halides,  is  far  from  extensive, 
it  should  disspell  the  common  view  that  simple  point  defects  are  the 
predominate  impurity  species  in  alkali  halides.  Furthermore,  most 
textbooks  create  the  impression  that  these  "simple  point  defects"  are  frozen 
in  some  reasonably  homogeneous  distribution  at  room  temperature.  There 
is  ample  evidence  that  this  is  just  not  true. 

At  high  temperatures,  an  alkali  halide  contains  most  of  its  impurities 

distributed  as  simple  point  defects,  but  their  geometric  distribution  is 

69 

not  uniform  enough  to  allow  a  measurement  of  the  Soret  effect. 

Even  if  quenched  to  low  temperatures,  dimers  and  trimers  of  impurity/ 
vacancy  complexes  are  very  quickly  the  predominate  impurity  species, 
and  in  three  months  to  a  year,  microprecipitates  will  have  formed 
and  OH"  and  H20  will  have  substantially  penetrated  the  crystal  if  it 
contains  a  few  parts  per  million  of  divalent  ions.  Annealed  crystals, 
even  if  relatively  pure,  contain  regions  of  abnormally  high  impurity 
concentration  because  clusters  and  precipitates  are  a  lower  total - 
energy  state  for  the  system. 

In  crystals  specially  purified  to  remove  impurity  cations  as  well 
as  anions,  the  highest  breakdown  thresholds  and  lowest  scattering  should 
be  observed  in  quenched  crystals,  provided  the  dislocation  density  can 
be  kept  reasonably  low.  However,  the  high  breakdown  thresholds  and  low 
scattering  would  be  reasonably  short  lived  because  smaller  impurity 
species  would  grow  slowly  as  the  crystal  ages. 

Molecular  ions,  with  their  permanent  dipole  moments,  are  expected 
to  be  strong  scatterers  because  they  are  strong  absorbers.  A  study 
of  molecular-ion  scattering  would  be  useful. 
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TABLE  III:  Experimental  Migration,  EM,  and  Aggration  E*  Energies 
and  Dimer  Formation  Energies  Eg  from  Crawfords  Model. 


System 

E*(ev) 

E„(eV) 

EB(eV) 

NaCl :Mn 

0.83  (30) 

0.95  (40) 

0.12 

0.78  ±  .05  (41) 

0.17 

NaCl :Cd 

0.50  ±  .04  (39) 

0.92  (41) 

0.42  • 

0.857  (43) 

0.36 

0.64  (42) 

0.14 

NaCl:Ca 

0.92  (41) 

0.96  (44) 

0.04 

0.90  (46) 

-0.02 

0.851  (45) 

-0.07 

1.15  (47) 

0.23 

NaCl :Pb 

0.80  (48) 

0.982  (43) 

0.18 

NaCl :Co 

0.61  (49) 

1.06  (50) 

0.45 

KC1 :Sr 

0.73  (30) 

0.871  (52) 

0.14 

KC1 :  Ba 

0.74  (30) 

KC1 :Pb 

0.85  (53) 

0.908  (54) 

0.06 
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TABLE  IV:  Enthalpies  and  Entropies  of  Solution  of  Divalent  Impurities 
in  Alkali  Halides. 


Host— - 
Impurity 


h(eV) 


KC1 

h(eV)  S(fc) 


h(eV) 


0.64 

-0.7(65) 

1.0 

(67) 

1.80 

15  (68) 

1.4 

(68) 

2.0 

9.4(68) 

1.6 

(68) 

0.086(27)' 


1.24  10.6(65)' 

1.14  (67)2 

0.992  0.81 (60) 1 

0.378  1.56(60)3,6 

1.195  0.95(60)5,6 

0.24  (27) (57)' 

0.816  -2.92(60)^ ,6 

0.386  -2.48(60)3,6 

0.347  -2.42(60)4,6 

0.984  -2.34(60)5,6 


1.14 

0.204 


0.190 


-  4.8(66) 


7.3(66) 


13  (66) 


Notes:  1.  Solute  assumed  to  form  M  and  V^  only  in  crystal. 

2.  hc,  the  enthalpy  of  formation  of  Vc  assumed  to  be  0.7  eV 

3.  Solute  assumed  to  form  ( MVC ( 0) } . 

4.  Solute  assumed  to  form  (MVC(1)}^. 

5.  Solute  assumed  to  form  ( MVC ( 0 ) } ^ . 

6.  Species  in  1,3, 4, 5  in  equilibrium  with  6NaCl  •  MCI.,. 
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D.  LIGHT  SCATTERING  BY  POINT  DEFECTS  IN  INSULATING  CRYSTALS* 


The  scattering  by  point  imperfections  and  clusters  is  often 
characterized  by  an  effective  particle  size,  which  is  obtained  formally 
by  equating  the  Mie-scattering  result  to  the  result  for  the  point 
imperfection  or  cluster.  This  effective  radius  rg  is  typically  greater 
than  the  actual  radius.  In  this  section  it  is  shown  that  the  scattering 
by  the  strain  field  around  an  impurity  ion  is  negligible  and  cannot 
account  for  the  difference  between  the  effective  radius  and  the  actual 
radius.  The  difference  between  the  effective  radius  and  the  actual 
radius  does  not  constitute  an  unresolved  theoretical  problem  because 
the  macroscopic  theoretical  Mie  results  is  not  expected  to  hold  at 
such  small  volumes  that  the  assumption  of  many  atoms  within  the  volume 
is  violated. 
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I.  INTRODUCTION 

There  is  currently  interest  in  the  elastic  scattering  of  light  from 

point  defects  in  alkali  halides,  since  these  materials  may  prove  useful 

for  fabrication  of  optical  fibers  with  ultra-low  absorption  near  1  pm^ . 

This  requires  that  all  sources  of  absorption  or  scattering  be  eliminated, 

save  that  from  intrinsic  processes.  In  materials  presently  available, 

quasi-elastic  scattering  of  light  is  observed  with  intensity  one  or  two 

orders  of  magnitude  greater  than  expected  from  theoretical  considerations 

applied  to  the  pure  material.  Experimentally,  one  verifies  this  by 

comparing  the  integrated  intensity  of  the  Brillouin  doublets,  and  that 

of  light  scattered  by  quasi -elastic  processes.  From  the  Landau-Placzek 

2 

ratio,  as  modified  recently  by  Wehner  and  Klein  ,  one  expects  the  quasi¬ 
elastic  peak  to  be  two  orders  of  magnitude  weaker  than  the  Brillouin 

doublets  in  alkali  halide  crystals,  while  in  practice  the  two  features  are 
found  to  have  comparable  intensity  even  in  crystals  of  high  nominal  purity. 

A  point  defect  in  the  crystal  scatters  light  elastically  because, 
when  its  influence  is  viewed  in  macroscopic  terms,  the  dielectric 
constant  within  the  unit  cell  that  contains  the  impurity  differs  from 
that  of  the  surrounding  medium.  This  leads  to  Rayleigh  scattering,  with 
cross  section  proportional  to  ai  ,  where  co  is  the  frequency  of  the  light. 
Thus,  point  defects  such  as  impurity  ions  or  vacancies  lead  to  scattering 
which  contributes  to  the  intensity  of  the  quasi-elastic  peak. 

In  addition  to  the  scattering  mechanism  described  above,  there  is  a 


second  means  by  which  even  an  isoelectronic  point  defect  may  scatter 
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light  elastically.  When  an  impurity  or  vacancy  is  introduced  into  the 
crystal,  the  nearest  neighbors  shift  off  lattice  sites  appropriate  to 
the  perfect  crystal,  with  the  result  that  a  macroscopic  strain  is  Induced 
in  the  surrounding  material.  This  strain  is  long-ranged,  and  has  the 
potential  of  substantially  increasing  the  cross  section  for  scattering 
from  the  defect.  When  it  is  noted  that  the  dielectric  constant  is 
affected  by  this  strain  field,  the  impurity  is  surrounded  by  a  "halo" 
of  substantial  effective  radius.  Even  if  the  local  change  in  dielectric 
constant  within  the  halo  is  small,  its  large  effective  radius  may  render 
it  effective  in  scattering  light. 

One  aim  of  this  paper  is  to  derive  an  expression  for  the  light 
scattering  cross  section  from  a  point  defect,  with  the  "strain  halo" 
effect  included.  Then  we  apply  the  resulting  formula  to  the  alkali 
halides,  with  emphasis  on  potassium  bromide,  which  has  been  studied 
extensively  by  Edwards  and  Fredericks3,  to  isolate  the  contribution  of 
particular  impurities  to  the  quasi-elastic  central  peak.  Our  simple  model 
is  applicable  to  isoelectronic  defects  in  alkai  halides,  and  here  we  reach 
conclusions  compatible  with  the  data:  that  isoelectronic  impurities  at  the 
part  per  million  level  fail  to  contribute  significantly  to  the  central 
peak. 


II.  THE  CALCULATION 


Let  the  defect  produce  a  perturbation  in  the  dielectric  tensor  of  the 
crystal  we  designate  as  Ae^Cx).  We  assume  that  the  crystal  is  optically 
isotropic,  but  the  strain  field  induces  local  anisotropy  in  the  dielectric 
tensor.  We  suppose  Ac.  .(x)  is  non-zero  only  in  the  neighborhood  of  the 

1  J 

defect,  although  we  may  proceed  without  specifying  its  form  in  detail. 

If  e|°^(x)  is  the  itfl  Cartesian  component  of  the  incident  electric 
field,  which  has  frequency  w,  then  to  first  order  in  Ac.  .(x)  the  scattered 

*  J 

field  is  given  by 


:js,(x)  =  -(f)2  f  d3x.£0lJ<;  -  x')Acjk(x’)Ek<°>(I') 
J  J  k 


(1) 


where  D.  .(x  -  x')  is  the  Green's  function  matrix  associated  with  Maxwell's 

'  J 

equations.  Far  from  the  impurity,  as  [x  |  with  x‘  fixed,  one  has 


D,.j(x  -  x-) 


.  ;<s>;<s>, -i ^(s)-x' 


4*  lx 


(5ij  - 


(2) 


1/2 

where  Kq  =  we  /c,  with  e  the  optical  frequency  dielectric  constant  of 
the  unstrained  crystal.  Here  we  assume  the  material  is  transparent,  so  e 
is  real  and  positive.  We  have  n'  '  a  unit  vector  from  the  impurity  to  the 
point  of  observation  x,  andt^  =  kQn^5^  is  the  wave  vector  of  the  scattered 
light.  If  for  e£°^(x')  we  take  the  form  6^  E^exp[i1c^0^*x],  corresponding 
to  an  electric  field  parallel  to  the  x-axis  with  wave  vector  1^°^  directed 

A 

parallel  to  z,  then  we  have 
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(x) 


g>2eiKox 

4ttc^x 


s«*n;  <•> 

®  V 


(3) 


where  we  have  defined 


Aei j (^)  -/A^' 


Aeij(x') 


(4) 


wi 


->  ■+($)  •+ 
th  Q  =  kv  '  -  k  the  scattering  vector.  The  magnitude  of  Q  i; 


2KQsin  (6/2)  where  0  is  the  angle  between  the  z-axis  and  the  wave  vector 
of  the  scattered  light. 

By  forming  the  Poynting  vector  associated  with  the  scattered 
radiation,  then  dividing  by  the  incident  energy  flux,  one  can  form  an 
expression  for  the  cross  section  a  for  scattering  from  the  point  defect. 
One  has 


(5) 


where  the  integration  is  over  all  directions  of  the  scattered  light. 

We  now  turn  to  our  explicit  model  of  Ac.  .(x).  Suppose  the 

*  J 

impurity  cell  perturbs  the  matrix  by  changing  the  dielectric  constant 
within  the  impurity  unit  cell  by  the  amount  Aeq.  If  V r  is  the 
volume  of  the  unit  cell,  approximated  by  a  sphere  of  radius  RQ,  we  then 
have 


=  +  4eu(*) 


(6) 
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where  Ac.:j(x)  is  the  part  that  comes  from  lattice  strains  outside  the 
unit  cell .  Thus 


-  Vc4c«(j  * 


.3  , 
d  x  e 


-i^*x 


(7) 


The  first  term  in  Eq.  (7)  taken  alone  gives  for  the  cross  section  the 
Rayleigh  result  oq  which,  in  our  present  model  is 


a 


o 


6irc^ 


(8) 


We  next  turn  to  the  form  of  Ae..(x).  We  model  the  strain  field  around 

*  J 

the  defect  as  follows.  Imagine  the  unit  cell  which  surrounds  the  defect 
is  a  sphere  of  nominal  radius  R  .  After  introduction  of  the  defect,  the 
impurity  cell  will  contract  or  expand,  to  change  the  radius  by  the  amount 
A.  This  produces  a  displacement  l5(x)  of  the  material  outside  the  unit 
cell  and  elasticity  theory  for  an  isotropic  material  allows  us  to  write 


where  x  is  the  unit  vector  in  the  radial  direction.  The  expression  in 
Eq.  (9),  and  of  course  the  results  deduced  from  it,  are  valid  only  in  the 
limit  where  Ais  small  compared  to  RQ,  so  linear  elasticity  theory  may  be 
used  to  describe  the  lattice  deformation  outside  the  unit  cell. 

In  the  presence  of  strains  described  by  the  strain  tensor 
Pij(x)  *  l/Zfai^./aXj  +  9Uj/3x.. ) ,  we  have  for  e^(x)  the  forms 
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■  - 


Pll' 


xx 


+  Pl2<eyy  +  ezz>  ' 


(10a) 


6£xyW  =  '  e2p44  exy  (10b) 

with  the  remaining  elements  of  6e..(x)  obtained  by  permutation.  In  Eqs.  (10), 
e  is  again  the  optical  frequency  dielectric  constant,  while  p^ ,  p^  and 
P44  are  the  Pockels  elasto-optic  tensors.  For  convenience,  we  adopt  the 
isotropy  assumption  p^4  =  p^  -  p^2  ,  which  allows  Eq.  (10a)  and 

Eq.  (10b)  to  be  combined  into  the  statements 

6£i/x>  ■  VijM)  ♦  -  V3  5(j(V.ti)|  (11) 

where  P]  =  -  e2(p11  +  2p12)/3  and  P?  =  e2(p^2  -  p^) 


Now  for  the  displacement  pattern  in  Eq.  (9),  one  has 

V*U  =  0  (12) 

and  there  is  no  local  dilation  in  volume  outside  the  defect  cell.  Thus, 

the  scattering  from  the  "strain  halo"  comes  entirely  from  the  second  term  in 

o 

Eq.  (11).  If  we  let  A  =  P2R£a,  ^en  91ves  us 


_3 

so  the  strain  field  falls  off  slowly  with  distance  as  x  .  With  this  form, 
one  has 
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4Mj(Q)  =c,(Q)SlJ  ♦e2(Q>Q1Qj 


(14) 


where 


4_A  /  sin  (QR0) 

ei  (Q )  =  V  Ae  + - 2  (  cos  (QR0} - 

1  c  0  (QR0)2  \  0  QR0 


(15a) 


and 


»(Q)  = 


127rft 

(QRq) 


/  sin  (QR  )  \ 

?(cos  <Q*V,I  — - ) 


(15b) 


For  our  purposes,  the  limit  QRq  «  1  is  appropriate,  so  we  may  replace 
these  forms  by 


E,(Q)  ■  Vco-T4 


(16a) 


and 


c2(Q)  =  +  4ttA 


(16b) 


It  is  a  straight  forward,  but  tedious  matter  to  use  Eq.  (14)  combined 
with  Eqs.  (16)  to  work  out  an  explicit  expression  for  the  light  scattering 
cross  section.  We  introduce  a  parameter  r  that  provides  a  dimensionless 
measure  of  the  role  played  by  the  "strain  halo".  We  thus  define 


r  = 


(17) 


which  enables  the  cross  section  to  be  written  in  the  form 
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°TOT  =  0oS 

where 


(18) 


S  =  1  -  |  r  +  |  r2  (19) 

and  where  the  quantity  aQ  is  the  Rayleigh  cross  section  defined  in  Eq.  (8). 

We  see  that  inclusion  of  scattering  from  the  strain  halo  alters  the 

4 

magnitude  of  the  Rayleigh  cross  section,  but  leaves  the  well-known  cj 
frequency  dependence  unaltered.  We  do  not  find  the  latter  result  obvious,  in 
view  of  the  long  ranged  nature  of  the  strain  field  around  the  defect.  In 
the  next  section,  we  estimate  the  value  of  r  for  a  typical  alkali  halide, 
along  with  the  scattering  power  of  a  defect  in  such  a  crystal. 

The  principle  results  of  the  present  paper  are  contained  in  Eq.  (18)  and 
Eq.  (19).  Since  these  results  were  obtained  within  the  context  of  a  simple 
model  which  treats  the  medium  outside  the  defect  as  an  isotropic  elastic 
continuum,  it  is  useful  to  comment  on  what  features  will  remain  in  a  more 
general  treatment  of  the  problem. 

It  can  be  shown  quite  generally,  even  if  the  material  outside  the  defect 
cell  is  not  istrotropic,  that  inclusion  of  the  scattering  from  the  "strain 

4 

halo"  leaves  the  oj  variation  of  the  Rayleigh  cross  section  unaltered.  This 
is  true  within  the  framework  of  a  calculation  that  treats  the  medium  outside 
the  defect  cell  as  an  elastic  continuum,  but  not  necessarily  isotropic  in 
character.  What  is  required  to  demonstrate  this  is  that  the  displacement 
field  t)( x )  falls  off  inversely  with  x2,  though  in  general  there  will  be  an 
angular  variation  of  the  displacement  pattern  more  complex  than  displayed 
in  Eq.  (9).  A  discrete  lattice  calculation  will  show  deviations  from  the 


Sec.  D 


79 


_2 

x  behavior  close  in  to  the  defect,  but  such  short  range  modifications  of  the 

4 

displacement  field  will  also  leave  the  a>  dependence  unaltered,  as  the  leading 
contribution  to  the  cross  section.  Thus  we  believe  that  Eq.  (18),  with  S 
independent  of  frequency,  to  be  valid  quite  generally  when  the  strain 
halo  scattering  is  incorporated  into  the  theory. 

Of  course,  the  specific  form  given  for  S  in  Eq.  (19)  is  dependent  on 
our  use  of  the  isotropic  continuum  picture  to  model  the  medium  outside  the 
defect  cell.  It  is  possible  to  obtain  formal  expressions  for  Ae-jj(tj)  for 
more  general  models,  and  in  specific  cases  it  should  be  possible  to  carry 
these  through  to  the  calculation  of  S.  But  these  forms  are  cumbersome  and 
not  enlightening,  so  we  have  confined  our  attention  to  the  simple  isotropic 
continuum  model,  which  we  believe  provides  reliable  semi -quantitative  estimates 
of  the  magnitude  of  S. 

We  have  also  confined  our  attention  to  an  isolated  defect  in  a  medium 
of  infinite  extent,  since  for  the  very  low  inpurity  concentrations  present 
in  the  high  purity  materials  that  motivate  this  work,  the  total  scattering 
cross  section  is  simply  the  cross  section  per  impurity  multiplied  by  the 
numbers  of  impurities  exposed  to  the  beam.  This  is  so  as  long  as  the 
impurities  are  distributed  randomly  in  the  host.  If  a  defect  is  placed 
in  a  finite  sample,  the  strain  field  must  be  synthesized  from  the  direct 
strain  field  from  the  defect  itself  (considered  here),  and  that  from  an 
array  of  image  defects  arranged  to  insure  satisfaction  of  the  stress  free 
boundary  conditions.  The  strain  fields  from  the  image  defects  are  small, 
and  will  lead  to  no  major  modifications  of  our  results,  save  for  defects 
located  near  the  surface.  In  highly  transparent  media,  these  are  only  a 
small  fraction  of  those  sampled  by  the  beam. 
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III.  DISCUSSION 


In  this  section,  we  comment  on  the  implications  of  the  result  displayed 
in  Eq.  (18),  which  is  the  principal  result  of  the  present  paper.  We  first 
assess  the  possible  importance  of  scattering  from  the  "strain  halo"  effect 
contained  in  the  factor  S  in  Eq.  (18). 

2 

For  a  number  of  crystals,  Wehner  and  Klein  have  tabulated  values 

of  p-j-j  and  p.j2,  along  with  the  optical  frequency  dielectric  constant.  This 
enables  us  to  calculate  the  constant  P2,  which  is  dimensionless  and  depends 
on  only  the  properties  of  the  host  crystal.  Both  A/Rq  and  AeQ  depend  on 
the  nature  of  the  particular  defect  under  consideration.  Here  we  choose 
|A/Ro|  =  0.1  and  |AeQ|  =  1.0  as  typical  order  of  magnitude  estimates  for 
these  quantities.  For  a  lattice  vacancy  in  an  alkali  halide,  where  the 
optical  frequency  dielectric  constant  is  near  2.0  in  value,  we  expect 
Agq  =  -  1.0  and  the  lattice  should  collapse  locally  into  the  defect  cell. 
Then  we  have  r  =  0.1  P2  for  this  picture. 

For  the  alkali  halides,  with  the  exception  of  lithium  flourine,  p-j , 
and  Dj2  do  not  differ  greatly.  For  potassium  bromide,  for  example, 

1  Pi  1  -  Pi 2 1  ~  0-04  so  with  c  =  2.36  we  estimate  | rj  ~  0.02  «  1.  The 
scattering  from  the  "strain  halo"  is  thus  quite  negligible. 

The  reason  for  this  is  clear  upon  noting  that  for  the  displacement 
pattern  in  Eq.  (9),  V*t5  =  0.  Thus,  outside  the  impurity  cell,  the 
lattice  distortion  leads  to  no  local  change  in  density,  but  it  simply 
lowers  the  local  symmetry  from  cubic  to  tetragonal.  The  dielectric 
constants  of  the  alkali  halides  is  described  rather  well  by  the  Clausus  - 
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Mcsotti  relation  combined  with  free  ion  polarizabilities;  this  is  a 
reflection  of  the  highly  ionic  character  of  the  alkali  halides.  The 
consequence  is  that  the  dielectric  constant  of  the  alkali  halides  depends 
on  the  density  of  the  ions,  but  is  insensitive  to  ionic  rearrangements 
which  leave  the  average  density  unchanged. 

The  above  argument  suggests  that  in  covalent  materials,  the  "strain 
halo"  scattering  may  be  much  more  significant  in  magnitude.  Consider 
silicon,  for  example.  Here  jp^  -  p^j  =  0,09,  and  P 2  =  12.9,  so  we 
estimate  |r|  -  1.29  with  the  above  value  of  (A/Rq)  and  Aeo.  Thus,  here 
the  "strain  halo"  scattering  is  estimated  to  be  of  equal  importance,  or 
possibly  larger  than,  that  produced  by  the  local  change  in  electronic 
polarizabil ity.  A  similar  estimate  for  barium  titinate  gives  | r|  =  1.42, 
so  again  the  "strain  halo"  scattering  contributes  importantly  to  the  cross 
section5. 

Our  conclusion  is  that  in  the  alkali  halides,  our  parameter  r  is 
small  compared  to  unity  a  1 wavs,  and  the  cross  section  provided  by  the 
Rayleigh  theory  is  quite  adequate.  In  less  ionic  materials,  such  as 
silicon  or  barium  titinate,  where  the  dielectric  tensor  is  sensitive  to 
local  changes  in  symmetry,  then  scattering  produced  by  the  strain  field 
outside  the  impurity  cell  is  of  comparable  importance  to  that  from  the 

local  change  in  electronic  polarizability  produced  by  the  defect. 

(3) 

Edwards  and  Fredericks'  have  carried  out  an  extensive  series  of 
light  scattering  studies  of  the  Landau- Placzek  ratio  in  potassium  bromide 
of  high  nominal  purity,  to  determine  whether  point  defects  are  responsible 
for  the  large  intensity  of  the  quasi-elastic  central  peak.  While  their 
samples  contain  rubidium  atthe  part  per  million  level,  present  as  an 
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isoelectronic  substitutional  impurity  which  replaces  potassium,  this 

impurity  has  no  effect  on  the  intensity  of  the  central  peak.  They  do 
2+ 

find  that  Fe  ,  present  at  much  lower  concentration  than  rubidium,  does 

2+ 

give  a  contribution  to  the  central  peak  linear  in  the  Fe  concentration. 

Our  model  applies  to  rubidium,  and  with  r  chosen  to  be  zero,  and 

Ae  =  -1,  we  find  the  Rayleigh  cross  section  assumes  the  value 

-25  2  °  + 

a  -  5.9  x  10  cm  .  This  estimate  applies  to  the  5145  A  Ar  laser 
o 

O 

line,  and  we  take  Rq  =  3.28  A,  the  nearest  neighbor  distance  in  potassium 

bromide.  If  rubidium  mpurities  are  present  at  the  one  part  per  million 

1  fi  3 

level,  then  with  n  -  1.4  x  10  defects  per  cm  ,  we  estimate  the  extinction 
length  associated  with  Rayleigh  scattering  from  rubidium  impurities  to  be 
£  =1.2  x  108  cm,  a  value  comfortably  higher  than  the  values  of  106  -  10^  cm 
typically  associated  with  Brillouin  scattering  intensities  in  the  alkali 
halides.  Thus,  our  conclusion  is  that  while  rubidium  is  the  most 

3 

abundant  impurity  in  the  samples  used  by  Edwards  and  Fredericks  ,  Rayleigh 
scattering  from  these  ions  is  too  weak  to  account  for  the  quasi-elastic  peaks 
observed  by  them.  This  is  in  accord  with  their  interpretation  of  the 
data,  as  outlined  above. 

2+ 

At  this  point,  we  can  say  little  about  scattering  produced  by  Fe 
on  the  basis  of  our  model,  which  ignores  the  macroscopic  electric  field 
outside  the  defect  unit  cell  present  when  an  impurity  with  valence 
different  than  the  host  ion  it  replaced  is  present.  Also,  such  impurities 
have  charge  compensation  centers  in  their  near  vicinity,  so  a  simple 
physical  picture  of  light  scattering  from  the  complex  involves  a  number 
of  considerations  not  contained  in  the  present  model. 


Sec.  D 


83 


IV.  SUMMARY 


We  have  examined  the  theory  of  light  scattering  by  point  defects  in 
crystals  through  use  of  a  model  which  includes  scattering  induced  by  the 
strain  field  near  the  defect,  in  addition  to  that  caused  by  the  change  in 
electronic  polarizability  within  the  unit  cell  which  contains  it.  The 
model  applies  to  defects  with  the  same  valence  as  the  missing  host  ion. 

We  find  the  light  scattering-cross  section  is  given  by  an  expression  of 
the  Rayleigh  form,  with  overall  magnitude  influenced  by  interference 
between  the  two  mechanisms  considered.  A  rough  estimate  of  the  cross 
section  shows,  in  agreement  with  recent  data,  that  isoelectronic  impurities 
present  in  alkali  halide  crystals  at  the  part  per  million  level  should  not 
influence  the  Landau-Placzek  ratio. 
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5.  In  silicon,  and  very  clearly  in  barium  titanate,  our  isotropy 
assumptions  must  be  modified  to  obtain  a  proper  description  of 
the  strain  halo  scattering. 


E.  MULTIPLE  SCATTERING 


In  this  section  it  is  shown  by  developing  a  multiple-scattering 
formalism  and  applying  the  result  to  scattering  by  randomly  located 
spherical  voids  that  the  multiple-scattering  reduction  of  the  individual 
scattering  results  is  entirely  negligible  in  low-loss  fibers  because 
the  scattering  must  be  small  and  multiple-scattering  effects  are 
negligible  for  small  scattering. 
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I.  INTRODUCTION 


The  propagation  of  light  in  a  medium  with  randomly  distributed 
inhomogeneties  or  scattering  centers  can  be  usefully  treated  by  the 
pertubation  methods  originally  developed  for  quantum  field  theory. 

In  addition  to  providing  a  convenient  framework  for  describing  the 
effects  of  multiple  scattering,  the  techniques  allow  important  classes 
of  terms  appearing  in  the  multiple  scattering  expansion  to  be  summed 
in  closed  form. 

Following  the  work  of  Tatarski\in  Sec.  II  we  derive  expressions 
for  the  ensemble-averaged  propagator,  which  carries  information  about 
the  effects  of  multiple  scattering  on  the  dielectric  properties  of  the 
medium.  In  particular,  it  provides  an  estimate  of  the  attenuation.  In 
Sec.  Ill  we  apply  the  results  to  the  simple  case  of  randomly  located 
identical  spherical  inclusions.  The  calculation  leads  to  simple  closed- 
form  expressions  for  the  modification  of  the  dielectric  constant  by 


the  inclusions. 
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II.  MULTIPLE  SCATTERING  FORMALISM 

Scalar  wave  propagation  in  a  spatially  varying  medium  is  de¬ 
scribed  by  the  Green's  function  G(r,  r')  between  points  r  and  r' 
satisfying  the  equation 

V2  G(r,  r’)  +  k2  e(r)  G(r,  r’)  =  6  3(r  -  r')  .  (1) 

Here  k^  5  2it/Aq  is  the  fr*_e  space  wave  number  corresponding  to  wave¬ 
length  Aq  and  c(r)  is  the  stochastic  dielectric  function  for  the  me¬ 
dium.  Defining  the  ensemble  averaged  dielectric  constant  <c>  and 
2  _  2 

^1  =  ^0  <'C'>  »  (1)  can  Ee  written 

V2  G(r ,  r’)  +  k2[l  +  E1(r)]  G(r,  r’)  =  63(r  -  r’)  (2) 

where 

e1(r)  H  (c(r)  -  <c>  ]/<t>  .  (3) 

2  2 

By  construction  <ep>=  0.  Defining  the  operators  L^  =  V  +  k^  and 
Mq  S  Lq  ^ ,  the  unperturbed  propagation  is  described  by  the  equation 
for  the  free  propagator  G^ 

L0  E’>  =  *3<E  -  £’) 


(4) 
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having  the  solution  satisfying  the  radiation  boundary  condition 
GQ(r,  r')  -  Mq  <S3(r  -  r’>  *=  -(47i)  1  exp  [ik^Jr  -  r'|]/|r  -  r'|.  (5) 

Substitution  into  Eq.  (2)  produces 

G(r,  r’)  =  -k*(MQ  C)(r,  r')  +  CQ(r ,  r')  .  (6) 

=  -k l  J d3r"  G0(r,  r")  e^r")  G(r",  r’)  +  GQ(r ,  r’).  (6a) 


Ensemble  averaging  Eq.  (6)  gives  in  operator  notation 

<G>  =  gq  -  MQ  <cx  G  >  .  (7) 

Equation  (7)  can  be  formally  solved  by  iteration  to  produce  the 
usual  Born  series.  Each  term  in  the  series  can  be  represented  by  a 
diagram  as  in  quantum  field  theory  leading  to  the  diagramatic  repre¬ 
sentation  of  Eq.  (7) 


+  .  .  .  + 
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Instead  of  solving  Eq.  (7)  by  iteration  directly  it  is  desirable 
to  separate  out  the  disconnected  diagrams  such  as 


Following  Tatarski  , 


set 


G  =  <  G  >  +  C„  (8) 

F 

so  <GI.>  =  0  and  substitute  into  Eq.  (2) 

r 

Lq[<G>  +  Gp]  +  kj  Ejl<G>  +  Gp]  =  63(r  -  r').  (9) 

Ensemble  averaging  gives 

LQ  <G>  +  k^  <  cx  GF  >  =  63 ( r  -  r')  (10) 


where  we  have  used  <C^>  =  0.  Subtracting  Eq .  (10)  from  (9)  and 
operating  with  Mg  produces 


CF  " 


-kj  Mq  Cj<G>  -  k^  MglCj  Gp  -  <c1  Gp  >  ] 


(ID 


Iteration  of  this  expression  followed  by  multiplication  by  k^  and 
ensemble  averaging  gives 


Unlike  Tatarski,  who  was  calculating  the  effects  of  normally  distri¬ 
buted  refractive  index  fluctuations  with  zero  mean,  we  cannot  in  gen¬ 
eral  neglect  the  odd  order  moments  of  e^, except  the  first.  This  leads 
to  the  appearance  of  terms  such  as  the  one  proportional  to  k^  which 
are  absent  in  Tatarski 's  treatment.  In  diagrams  the  expression  above 
for  Q  appears  as 


This  represents  the  "self-energy"  contribution  to  the  propagator  and 
consists  of  only  the  connected  diagrams  of  <G>  .  Substituting 
Eq.  (12)  into  (10)  gives 

LQ  <G>  -  Q  <  G  >  =  63(r  -  r’).  (13) 

Operating  with  produces 

<  G  >  =C0  +  M0Q<C>  (14) 
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which  is  Dyson's  equation  representing  the  average  propagation  as  an 
infinite  sum  of  "self-energy"  terms  coupled  by  the  free  propagation. 
In  diagrams 


■ - +  -O— +  —0—0— *7 

Equation  (14)  can  be  solved  using  Fourier  transforms  once  a  choice 
of  Q  is  made;  letting  <g>  ,  gQ ,  and  q  devote  the  respective  trans¬ 
forms  of  <  G  >  ,  Gq ,  Q  a  simple  calculation  gives 


<  g(f)  >  =  gQ(K)  +  (271)6  q(K)  <g(<)>  (15) 

or 

<g(<)>  =  (g0(K)_1  -  (  2tt  ) 6  q  (<)  ]  ~ 1  .  (16) 


If  desired,  <G>  can  be  obtained  from  Eq .  (16)  by  inverse  transforma¬ 
tion;  however,  the  most  useful  information  about  the  averaged  propaga¬ 
tor  concerns  the  modification  of  multiple  scattering  to  the  dielectric 
constant  which  can  be  obtained  directly  from  the  pole  structure  of 
Eq.  (16)  as  will  be  illustrated  in  the  next  section. 
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III.  SCATTERING  FROM  SPHERICAL  INCLUSIONS 


As  an  example  of  the  application  of  the  formalism  just  developed, 
consider  the  scattering  from  a  randomly  distributed  set  of  identical 
spherical  inclusions  of  radius  a  and  dielectric  constant  imbedded  in 
a  uniform  background  of  dielectric  constant  £q.  The  dielectric  function 
for  a  sample  volume  ft  is  then 

N 

e(r)  W(r  -  R.)  (17) 

i  =1 

where  we  have  assumed  there  are  N  inclusions  located  at  random  points 
R ^ .  In  this  case 


W(r)  = 


for  |  r  -  R..  |  <  a 


for  |r  -  R^|  >  a  (18) 


The  mean  of  e(r)  is 


N  N 

<  e  >  =  ff  f  =  sT  2  /  W(r  =  R i)  =  2  /  d3p  w(^} 

i  =  1  i  =  1 


-  fes  +  (1  -  f)eQ 


(19) 


where  f  =  4ira  N/3ft  is  the  fractional  volume  displaced  by  all  the  inclusions. 


A 
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To  estimate  Q  we  make  the  simplest  approximation  of  keeping  only  the 
lowest  order  contribution  and  assume  the  spheres  do  not  overlap.  Then 

Q(r  -  r’)  ^  <  Cj  MQ 

-  kj  £j(r')>  G0(r  -  r') 

=  - 5-  t  <£(r)  e(r')>  -  <c  >  ]  G  (r  -  r').  (20) 

<e>  u  ” 


Let  Be(r  -  r ' )  =  <  e(r)  e (_r ' )  >  be  the  correlation  function  satisfying 
the  requirements  Be(0)  =  8c j:  +  (1  -  B)  £g  and  r^moo  B£(_r)  =  (e  )2 
Expanding  <  e  >,  it  can  be  shown  that 

B£(r)  =  6(1  -  6)(es  -  e0)2  f(r)  +  <c>  2  ,  (21) 

where  f(r)  is  the  overlap  function  for  two  spheres  of  radius  a.  Straight¬ 
forward  calculation  gives 


f(r) 


,  r  <  2a 


0 


r  >  2a  . 


(22) 


«  j*  /  a 

The  mathematics  is  greatly  simplified  by  approximating  f(r)  =  e  for 
all  r.  With  this  choice 
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B£  (r)  ~  3(1  -  6)(cs  -  eQ)2  e"r/a/<c>2 


"~3  B 

!tt)  J  J  1 


(r)  G0(r)  e  ^ 


ko  (cs  "  co)2 

- 3  6(1  -  B)  --  -■  —  - 2 

(2tt)  k  +  (-  -  ik.r 

H  1. 


Substituting  into  Eq.  (16)  gives 


<g(<)>=  (  2tt  ) 


k2  +  k*  3(1  -  6)(es  -  e/ 
k2  +  (a  "  ikl)2 


For  small  particles,  with  k^a  <<  1,  the  poles  of  <gU)>  occur  at 


k2  -  k2  +  kj  3(1  -  B)(es  -  e0)2  a2  (1  +  i2k  a)  . 


Defining  k-j  =  pQ  and  k  =  ngp0  +  iBsc>  we  find 


ncc  [l  +  k2  a2  3(1  -  3)  Re(rs  -  LqY 2/n^  >  ] 1/2 
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3 


sc 


kg  a3  f(l  -  f)  Re(e$  -  Cg)2 


(28) 


Here  ng  represents  the  correction  to  the  average  refractive  index  n^  £  >  , 
defined  as  Pg  =  k^n^  £  y,  due  to  multiple  scattering  from  the  inclusions 
and  $sc  is  the  scattering  coefficient  for  Rayleigh  scattering  corrected 
for  multiple  scattering  through  the  factor  n  . 

The  numerator  of  the  second  term  is  essentially  the  Rayleigh 

2 

scattering  result,  given  by  Condon  and  Odishaw  for  example.  Since 
ne  >  1,  multiple  scattering  has  the  effect  of  decreasing  the  single 
scattering  contribution  to  the  extinction  coefficient  through  rescat¬ 
tering  into  the  forward  direction. 

To  estimate  the  importance  of  this  screening  we  assume  f  «  1  and 
approximate  Eq.  (27)  as 


n  =  1  + 


(k0ar 


Re(e, 


-'0>‘ 


(29) 


Tn  order  for  the  Rayleigh  model  to  be  valid,  the  inequality  kQa  £  tt  must 
be  satisfied.  Consider  the  extreme  case  of  kga  =  n  and  scattering  by 
vacuum  inclusions  (e^  =1).  Figure  1  is  a  plot  of  np  as  a  function  of 
the  refractive  index  ng  of  the  host  material,  taken  to  be  real,  for 
various  values  of  f.  We  see  that  in  order  to  observe  a  ten  percent 
reduction  in  the  extinction  coefficient  for  ng  =  2,  f  must  be  about 
10  ;  since  f  can  be  expressed  as  (a/s)  where  s  is  the  average 

spacing  between  inclusions,  this  implies  a/s  =  0.2  which  is  quite  large. 
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We  conclude  that  for  a/s  «  0.2  multiple  scattering  is  not  likely  to 
change  the  first  order  Rayleigh  result  appreciably.  For  the  case  of 
low-loss  fibers,  the  multiple-scattering  reduction  of  the  individual 
scattering  results  is  entirely  negligible  because  the  scattering  must 


FIGURE  CAPTIONS  FOR  SEC.  E 


Refractive -Index  dependence  of  correction  factor  ne  for 
multiple  Rayleigh  scattering  from  spherical  vacuum  inclusions 
imbedded  in  a  material  of  refractive  index  nQ.  The  parameter 
f  represents  the  fractional  volume  displaced  by  the  inclusions. 


Average-refractive-index  correction,  nfi  (dimensionless) 


Sec.  E 


99 


Unperturbed  refractive  index,  Op  (dimensionless) 


Fig.  1— Refractive-index  dependence  of  correction  factor  ne 

for  multiple  Rayleigh  scattering  from  spherical  vacuum 
inclusions  Imbedded  in  a  material  of  refractive  index 
no-  The  parameter  f  represents  the  fractional  volume 
displaced  by  the  inclusions. 
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F.  TABULATION  OF  INFRARED-ABSORPTION  COEFFICIENTS 


The  lowest  experimental  values  of  infrared  optical-absorption 
coefficients  from  the  literature  are  tabulated  as  plots  of  absorption 
coefficients  as  functions  of  frequency.  Beth  intrinsic  and  extrinsic 
values  are  tabulated. 


ft  P 1 y  1  °  0 ;  j  7 
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There  has  been  considerable  recent  interest  in  lowering  the  ab- 
sorptance  of  transparent  materials  for  use  in  high-power  laser  appli¬ 
cations.^  Both  the  intrinsic  values  of  the  optical-absorption  coef¬ 
ficients  and  the  lowest  experimental  values  are  useful.  Early  in  the 

2 

study  of  high-power  infrared  laser  materials.  Sparks  tabulated  absorp¬ 
tion  coefficients  that  were  available  in  the  literature  at  that  time. 
These  absorption  coefficients  were  extremely  useful,  but  that  tabulation 
now  includes  only  a  fraction  of  the  measured  absorption  coefficients 
because  most  of  the  currently  available  data  have  been  taken  since  the 
original  tabulation  was  made.  Other  properties  of  laser  materials  that 

have  been  tabulated  include  the  following:  ultraviolet  and  visible- 

3 

impurity  absorption  spectra,  extrinsic  absorptance  in  infrared  laser- 

4  5  6 

window  materials,  ’  and  laser-damage  thresholds. 

Here  we  tabulate  the  lowest  experimental  values  of  the  optical 
absorption  coefficients  3  reported  to  date  for  infrared  transparent 
materials.  The  results  for  all  materials  are  plotted  on  the  same  scale 
for  convenience  of  direct  comparison.  Since  the  multiphonon  absorption 
mechanism  is  now  well  understood,^  the  curves  are  more  useful  in 
practical  applications  than  in  theoretical  investigations.  Thus,  the 
absorption  coefficients  are  plotted  as  functions  of  the  frequency  w, 
rather  than  the  ratio  w/tOf,  where  uy  is  the  fundamental  Reststrahlen 
frequency.  The  strong  frequency  dependence  of  the  real  part  of  the 
index  of  refraction  near  the  fundamental  resonance  at  w  =  Wf  causes 
some  inaccuracy  of  the  values  near  the  fundamental  Reststrahlen 
resonance  peaks. 

_2  - 1 

The  larger  values  of  6,  typically  greater  than  10  cm  ,  can  be 
measured  spectroscopically.  Measurement  of  smaller  values  of  6 


requires  calorimetric  measurements,  for  which  only  a  limited  number 
of  discrete  laser  frequencies  have  been  available  in  the  measurements 
to  date.  Other  techniques  of  measuring  small  absorption  coefficients, 
such  as  emissivity,  frequency-modulation  spectroscopy,  and  photo- 
acoustical  measurements,  have  been  demonstrated  but  not  used  extensively. 

The  general  shape  of  the  3{<d)  absorption  curves  is  shown  sche¬ 
matically  in  Fig.  1.  The  values  of  3  on  the  frequency  range  from  the 

-2  -5  -1 

Reststrahlen  absorption  peaks  to  the  values  of  3  as  10  to  10  cm 
(depending  on  the  material  and  its  purity)  are  generally  believed  to  be 
intrinsic--that  is,  a  property  of  a  perfect  sample.  For  example,  the 
absorption  coefficient  for  potassium  chloride  in  Fig.  22  for  fre¬ 
quencies  to  <  1000  cm"l  is  believed  to  be  intrinsic.  The  value  of  the 
intrinsic  absorption  coefficient  decreases  as  the  frequency  increases 
for  id  >  iDf.  This  strong,  nearly  exponential  frequency  dependence  is 
well  explained  by  multiphonon  absorption.7  The  intrinsic  multiphonon 
absorption  becomes  so  small  at  high  frequencies  that  extrinsic  processes 
dominate  the  absorptance.  For  example,  at  frequencies  to  >  1000  cm~^ 
for  potassium  chloride,  the  absorptance  is  extrinsic  (up  to  the  in¬ 
trinsic  electronic-absorption  edge). 

4 

The  source  of  the  extrinsic  absorption  is  seldom  known.  Duthler 

has  shown  that  molecular-ion  impurities  are  a  major  potential  source 

8  9 

of  impurity  absorption.  Hooper  and  Uhlmann  and  Sparks  and  Duthler 

showed  that  macroscopic  inclusions  were  important  sources  of  intrinsic 

absorption. 

The  conductivity  a,  the  extinction  coefficient  k,  and  the  imaginary 
parts  Cj  of  the  dielectric  constant  are  given  in  terms  of  the  absorp¬ 
tion  coefficient  8  as  follows: 
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a  =  (cnr/4TT)B  , 
k  =  (X/47r)g  , 

£j  =  (cnr/o)}8  , 

where  u>  is  the  laser  frequency,  c  is  the  velocity  of  light,  and  nr  and 
the  extinction  coefficient  k  are  the  real  and  imaginary  parts  of  the 
index  of  refraction,  n  =  nr  +  ik. 

The  data  in  Figs.  2  through  53  are  from  references  10  through  64. 
The  curves  and  data  points  are  marked  with  the  reference  numbers. 
Materials  in  Figs.  2  through  53  are  listed  in  Table  I. 

We  greatly  appreciate  the  assistance  of  Mrs.  Lona  Case  and  Mr. 

D.  F.  King  in  preparing  the  tables.  This  research  was  supported  by 
the  Defense  Advanced  Research  Project  Agency,  administered  by  the 
Naval  Research  Laboratory,  the  Air  Force  Office  of  Scientific  Research, 
and  the  Air  Force  Wright  Aeronautical  Laboratories. 
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TABLE  I 


List  of  materials  included  in  the  tabulations. 


AgC^ 

KCS- 

m2o3 

KI 

MSb 

KMgF3 

^s2^3 

KRS-5 

As2Se3 

KRS-6 

BaF2 

KZnF3 

BaTi03 

LaF2 

BeF2 

Li  F 

CaC03  (calcite) 

LiYF4 

CaF2 

MgF2 

CdTe 

MgO 

CsBr 

NaQ 

Csl 

NaF 

Diamond 

Nal 

GaAs 

RbU 

Gd3Ga51012 

Si 

(gadolinium-gallium 

garnet) 

bl  L 

Ge 

Si 02  (crystalline 

quartz) 

HfF*  (based  ternary 

glass) 

Si02  (fused  silica) 

InAs 

SrF2 

KBr 

SrTi03 
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Fig.  1. 


Figs.  2 


FIGURE  CAPTIONS  FOR  SEC.  F 


Schematic  illustration  of  the  frequency  dependence  of  the 
absorption  of  a  dielectric  material,  showing  the  intrinsic 
multiphonon-  and  electronic-absorption  and  the  extrinsic 
absorption  range. 

-  53. 

Lowest  experimental  values  of  absorption  coefficients  at 
various  frequencies. 
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G.  V  CURVES  OF  ABSORPTION  AND  SCATTERING 


Plots  of  the  intrinsic-absorption  coefficients  and  the  intrinsic- 

scattering  coefficients  have  been  called  V  curves.  Such  curves  are 

useful  in  displaying  the  intrinsic  limits  of  the  loss  in  fiber.  However, 

they  are  likely  to  be  misleading  because  neither  the  intrinsic  scattering 

limit  nor  the  intrinsic  absorption  limit  is  currently  attainable  in  most 

cases  of  current  interest.  The  V  curves  are  more  useful  in  establishing 

the  frequency  region  in  which  the  intrinsic  absorption  and  scattering  will 

allow  operation  at  a  specified  loss  (8gx  =  Bab  +  Bsc)  than  in  determining 

the  minimum  of  3ex  and  the  corresponding  wavelength.  For  example,  if  an 

-8  -1 

extinction  coefficient  Bgx  =  +  Bsc  =  10  cm"  is  required,  the 

-8  -1 

useful  region— in  which  8  10”  cm  —for  KI  is  approximately  1,900  cm" 

CA 

<  a)  <  4,800  cm”^  (or  2.1  ym  <  X  <  5.3  ym) ,  as  seen  in  Fig.  8.  It  is  not 
necessary  to  operate  at  the  minimum  of  the  B  (o>)  curve,  which  occurs 

Ca 

near  2,300  cm~^  (or  4.4  ym)  for  KI,  as  seen  in  Fig.  8.  Inaccuracies  in 
extrapolating  the  absorption  curves  and  in  estimating  the  intrinsic 
scattering  make  the  operating-region  limits  quite  approximate.  With  these 
important  reservations,  the  V  curves  in  Figs.  1  through  8  are  presented. 
The  absorption  curves,  which  establish  the  low-frequency  side  of  the  V, 
are  obtained  by  extrapolating  the  intrinsic  portions  of  the  Bab(w)  curves 
of  Sec.  F  to  high  frequencies.  The  solid  portions  of  the  absorption 
curves  are  in  the  region  in  which  the  measurements  were  made,  and  the 
dashed  portions  are  extrapolations.  A  few  of  the  absorptance  curves 
are  even  more  unreliable  because  they  are  extrapolations  (dotted)  of 
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measured  curves  (dashed-dotted)  that  are  too  short  to  extrapolate  with 
confidence.  The  Brillouin-scattering  curves,  which  are  all  plots  of 
theoretical  results,  are  solid  if  all  the  parameters  in  Eq.  (1)  below 
are  known  and  are  dashed  if  any  parameters  (such  as  p^)  are  obtained 
by  formally  using  values  for  similar  materials. 

The  Brillouin-scattering  curves  for  the  crystalline  materials  are 
plots  of  =  B$c,  with 

BB  =  8Ar8p122kBTB/3X4  (1) 

where  p^2  is  a  strain-optic  coefficient,  kg  is  the  Boltzman  coefficient, 

1  ? 

and  B  is  the  isothermal  compressibility.  ’This  result  is  for  both  the 
incident-  and  scattered-radiation  polarized  along  the  z  axis.  In 
general,  p^B  is  replaced  by  different  combinations  of  the  strain-optic 
coefficients  and  elastic  constants,  depending  on  the  polarizations  and 
measurement  geometry.  The  order  of  magnitude  of  the  resulting  value 
of  Bg  does  not  change  in  general. 

The  intrinsic  scattering  curves  for  glasses  are  obtained  as  follows: 

2 

The  curves  for  Si02,  ZnC£2,  and  BeF2  are  from  van  Uitert  and  Wemple. 

The  curves  for  As2Sg  and  TI  1173  glass  are  obtained  from  Eq.  (1)  with  T 
replaced  by  an  effective  temperature  which  is  approximated  by  the  glass 
softening  temperature,  573  K  for  As2S3  and  643  K  for  TI  1173  glass. 

The  corresponding  values  for  SiO^  and  ZnCf2  are  1700  K  and  556  K,  respec¬ 
tively.  Formally  replacing  the  temperature  by  the  effective  temperature 
accounts  for  the  random  molecular  structure  of  the  glass,  which  is 
visualized  as  phonon-type  fluctuations  that  are  frozen  into  the  structure 
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when  the  glass  solidifies.  As  an  example  of  the  magnitude  of  the  random- 

molecular-structure  scattering,  the  scattering  coefficient  for  SiC^  is 

a  factor  of  1700/293  =  5.8  times  greater  than  the  Brillouin  scattering. 

The  zero-dispersion  crossover  wavelengths  (at  which  d  n/dA.  =  0) 

are  marked  on  the  figures.  The  values  were  obtained  from  the  minima  of 

2  2 

the  curves  of  dn/dA  in  references  3  and  4.  Notice  that  d  n/dw  is  not 

2  2 

zero  at  the  same  wavelength  at  which  d  n/dX  =  0  in  general. 

Values  of  the  parameters  used  in  the  calculation  are  listed  in  Table  I. 
The  frequency  dependence  of  the  stress-optic  coefficients  p^  are  neglected. 
This  is  a  very  good  approximation  because  the  scattering  is  important  at 
short  wavelengths,  X  <!  5  pm  typically,  and  the  wavelength  dependence  of 
p^2  at  short  wavelengths  is  quite  weak. 
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Table  I.  Values  of  parameters  used  in  calculating  from  Eq.  (1). 


Material 

n 

(dimen¬ 

sionless) 

p12 

(dimen¬ 

sionless) 

V 

(dimen¬ 

sionless) 

E 

(106psi) 

(5) 

bt 

& 

Bg(2.8  pm) 

(10'10  cm"1 

SiO^ 

1.45 

9 

0.281 

a 

0.376  k (3] 

n.o  i 

0.981 

8.23 

^s2^3 

2.41 

b 

0.272 

e 

0.3  (2) 

2.3  k 

7.575 

3,470 

TI  1173 

2.60 

c 

0.213 

a(i: 

0.25  k 

3.1  l 

7.023 

3,620 

A^2°3 

1.67 

c 

0.038 

e 

0.248  k (3) 

50  l 

0.439 

0.209 

BaTiOg 

2.4 

f 

0.1 

i 

0.267  n 

23.5  n 

0.863 

51.7 

MgO 

1.67 

c 

-0.04 

e 

0.232  k (3) 

36.1  l 

0.646 

0.340 

Diamond 

2.42 

b 

0.325 

e 

0.29  g(3) 

100  l 

0.183 

123.7 

Ge 

4.02 

b 

-0.130 

a 

0.273  k(3) 

14.9  l 

1.326 

8,187 

Si 

3.43 

c 

0.01 

e 

0.278  k(3) 

19.0  l 

1.017 

10.6 

CdTe 

2.69 

c 

-0.017 

d 

0.408  k(3) 

5.3  L 

1.512 

1.92 

GaAs 

3.30 

c 

-0.140 

e 

0.311  k(3) 

12.3  l 

1.344 

594 

ZnS 

2.25 

a 

-0.01 

d 

0.385  k(3) 

7.73k(4) 

1,216 

0.128 

ZnSe 

2.43 

a 

-0.04 

a 

0.376  k(3) 

10.3  t 

1.048 

3.27 

BaF2 

1.45 

c 

0.194 

a(l) 

0.31  k(3) 

7.7  l 

2.148 

8.59 

CaF  2 

1.41 

a 

0.095 

e 

0.23  k(3) 

11.0  l 

2.133 

1.64 

MgF  2 

1.38 

c 

0.12 

(2) 

0.382  k(3) 

24.5  l 

0.419 

0.43 

SrF2 

1.42 

a 

0.143  a(l) 

0.258  k(3) 

14.7  l 

1.434 

2.63 

NaF 

1.32 

b 

0.20 

d 

0.202  k(3) 

1 2 . 6  k  (4 ) 

2.060 

4.13 

LiF 

1.39 

b 

0.130 

e 

0.272  g(3) 

9.3  k 

2.133 

2.70 

Nate 

1.52 

c 

0.165 

e 

0.252  j 

5.8  l 

3.729 

15.8 

KC€ 

1.47 

a 

0.15 

a 

0.146  k ( 3 ) 

4.3  t 

7.169 

19.1 

KBr 

1.52 

b 

0.171 

e 

0.144  k(3) 

3.9  l 

7.952 

36.0 

KI 

1.67 

b 

0.169 

e 

0.138  g(3) 

3.7  m 

8.518 

80.0 
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Notes: 

(1)  Value  of  P.| 2  calculated  from  and  q12. 

(2)  Value  estimated  from  related  materials. 

(3)  Value  of  v  calculated  from  v  =  C-j^/CC-^  +  C12). 

(4)  Value  of  E  calculated  from  E  =  (C^  +  2C12)(C-j1  -  C-|2^C11  +  C12). 

(5)  Value  of  By  calculated  from  =  3(1  -  2v)/E. 
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FIGURE  CAPTIONS  FOR  SEC.  G 


Fig.  1.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic, 
random-molecular-structure  scattering  coefficient  for  the  glasses  ZnCf2, 
SiO^,  and  BeF2.  In  this-  and  the  following-figures,  the  solid  portions 
of  the  absorption  curves  are  in  the  region  in  which  the  measurements  were 
made,  and  the  dashed  portions  are  extrapolations.  Dotted  absorption  curves 
are  extrapolations  of  measured  curves  (dash-dotted)  that  are  too  short  to 
extrapolate  with  confidence.  The  theoretical  Brillouin-scattering  curves 
are  solid  if  all  parameters  are  known  and  are  dashed  if  any  parameters 
are  obtained  by  formally  using  values  for  similar  materials. 

Fig.  2.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic, 
random-molecular-structure  scattering  coefficient  for  the  glasses  As2S3 
and  TI  1173.  See  the  caption  of  Fig.  1  for  the  meaning  of  the  various 
1 i nes . 

Fig.  3.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  A£203,  MgO,  and  BaTi03- 
See  the  caption  of  Fig.  1  for  the  meanings  of  the  various  lines. 

Fig.  4.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  Ge  and  Si.  See  the 
caption  of  Fig.  1  for  the  meanings  of  the  various  lines. 

Fig.  5.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  CdTe,  GaAs,  ZnS,  and 
ZnSe.  See  the  caption  of  Fig.  1  for  the  meanings  of  the  various  lines. 
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Fig.  6.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  BaF^,  CaF2,  MgF2  and 
SrF2.  ^ee  caPt’’on  °f  Fig.  i  for  the  meanings  of  the  various  lines. 

Fig.  7.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  NaF,  Li F ,  and  NaC£. 

See  the  caption  of  Fig.  1  for  the  meanings  of  the  various  lines. 

Fig.  8.  Extrapolated  intrinsic  absorption  coefficient  and  intrinsic 
Brillouin  scattering  coefficient  for  the  crystals  KC l,  KBr,  and  K1 .  See 
the  caption  of  Fig.  1  fo  the  meanings  of  the  various  lines. 
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H.  EXTRINSIC  ABSORPTION 


In  this  section  it  is  shown  that  only  very  small  trace  amounts  of 
either  absorbing  impurities  {as  low  as  10-4  parts  per  billion)  or  macro- 
scopic  inclusion  (as  low  as  10  parts  per  billion)  can  give  rise  to  the 
typical  goal  value  of  3ex  =  10  dB/km.  Tabulations  of  molecular-ion 
absorption  information  given  in  Table  I  and  Figs.  1  through  4  and  of 
visible-  and  ultraviolet-absorption  spectra  should  be  useful  in  identifying 
troublesome  impurities. 

The  concentration  of  absorbing  impurities  must  be  kept  at  extremely  low 

values  in  order  to  attain  even  the  typical  goal  value  of  Bgx  =  10  dB/km. 

In  recent  high-power-window  investigations,  research  programs  were  required  for 

most  materials  to  reduce  the  absorption  coefficient  to  Bgb  =  10'4  cm'1. 

Caloremetric-measurement  techniques  had  to  be  developed  in  order  to 

_2  -1 

measure  values  of  6ab  1  ess  than  approximately  10  cm  .  Apart  from 

reducing  the  absorptance,  even  the  measurements  of  values  of  the  absorption 

-4  -1 

coefficient  less  than  10  cm  has  been  difficult.  These  results  put  the 

-2  -8  -1 

typical  goal  value  of  8ab  =  10  dB/km  =  2.3  *  10”  cm”  in  perspective. 
Molecular-ions  are  believed  to  be  a  major  source  of  extrinsic  infrared 

'I  O 

absorption,  as  shown  by  Duthler.  Deutch  showed  that  molecules  on  surfaces 

of  windows  gave  rise  to  absorption  at  the  level  of  absorptance  A  =  10'3  to 
-4 

10  .  Spectroscopists  have  measured  the  absorption  spectra  of  molecules  for 

years  by  doping  then  into  alkali  halide  crystals.1  These  studies  provide 

a  wealth  of  information  on  absorption  by  the  molecular-ion  impurities. 

1  3 

Duthler  and  Flannery  and  Sparks  tabulated  such  information  that  related  to 
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impurity  absorption.  The  results  of  Flannery  and  Sparks'3  are  repeated 

here  in  Figs.  1  through  4  and  Table  I  for  the  convenience  of  the  user. 

Nominally  pure  commercial  crystals  are  commonly  observed  to  contain 

-  4  5 

molecular-anion  impurities,  such  as  the  metaborate  ion  (BO?)  ’  and  the 
2-  3 

carbonate  ion  (CO^  )  ,  using  conventional  absorption  spectroscopy.  The 
hydroxyl  ion  (0H~)  is  especially  difficult  to  remove  from  alkali  halides 
during  purification  ,  and  may  be  the  ion  limiting  absorption  at  2.7  Pm. 
These  anion  impurities  along  with  divalent  metal  cation  impurities6,7  are 
most  likely  still  present  in  small  concentrations  in  specially  purified 
materials  because  concentration  at  the  parts  per  billion  level  can  give 


rise  to  the  presently  common  extrinsic  values  of 


=  10"4  cm  ^ . 


The  severity  of  the  problem  of  the  limit  on  the  loss  3gx  by  trace 
amounts  of  molecular  ions  can  be  seen  from  the  results  of  Duthler1  and  of 

3 

Flannery  and  Sparks.  Scaling  the  result  that  parts  per  billion  can  give  rise  to 
6  .  =  ICT4  cm-1  to  the  typical  goal  of  e  =  1 0~2  dB/km  =  2.3  *  10~8  cm-1 

3D  6X 

-4 

shows  that  10  parts  per  billion  of  absorbing  impurities  can  give  rise  to 
10"2  dB/km. 

This  result  can  also  be  seen  directly  from  the  following  general 

5  -1 

argument.  From  the  typical  value  of  £ab  =  10  cm  at  Reststrahlen 

22 

absorption  peaks  in  solids,  with  typical  densities  of  Nsqj  =  10  ions 
3 

per  cm  ,  a  typical  infrared  absorption  cross  section  of  a  strong  absorber 
at  resonance  is 

°ab  =  6ab/Nsol  =  10  cm  710  cm  =  10'  cm\  (1) 

_2 

The  corresponding  concentration  N  required  to  give  rise  to  Bgb  =  10  dB/km 
=  2  x  10-8  cm  ^  is 

N  =  3ab/ogb  =  2  x  10~8  cm-1/10~17  cm2 

=  2  x  109  cm'3  (2) 
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22  3  9-3 

For  ~10  ions/cm  in  solids,  the  impurity  density  N  =  2  *  10  cm 

corresponds  to  the  concentration  2  *  10^/10^  =  2  x  10^  =  2  x  10"^  parts 

per  billion.  Even  the  measurement  of  such  low  concentrations  of  important 

impurities  is  impossible  at  present. 

8  9 

>  Another  source  of  extrinsic  absorption  is  macroscopic  inclusions.  ’ 

Platinum  particles  in  laser  glass  is  a  well  known  example.  Other  examples 

include  nonstoichiometric  regions  in  crystals,  contaminated  grain  boundaries 

g 

or  dislocations,  and  clusters  of  impurities.  Sparks  and  Duthler  showed 

-8 

that  fractional  volumes  as  low  as  10  can  give  rise  to  an  absorption 
coefficient  of  8ab  =  10”^  cm""*.  Formally  scaling  linearly  gives  a  fractional 
volume  of  2  x  lo~12  =  2  x  10~8  parts  per  billion  for  Bab  =  10  ^  dB/km 
=2x  10"8  cm"1. 
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Figure  2.  Correlation  chart  showing  the  absorption  lines  between  0.5  and 
15  pm  of  some  chemical  bonds,  molecular  Ions,  radicals,  and  compounds. 

•  peak  positions 

-  range  of  peak  positions  In  various  materials 

Absorption  strength: 

S  strong 
M  medium 
W  weak 
V  variable 
Absorption  line  width: 

B  broad 
Sh  sharp 

Numbers  associated  with  an  absorption  band  indicate  the  number  of 
peaks  occurring  within  that  band. 

If  several  designations  apply  to  one  peak,  they  are  separated  by 
conmas. 

If  a  designation  applies  to  two  adjacent  peaks,  it  may  be  written 
between  them. 
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Alkanes,  cycloalkanes 

Alkenes,  cycloalkenes 

Alkynes 

Aroaatlcs 

Fluorides 

Chlorides 

Broaldes 
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Alcohols 

Phenol s 

Ethers 

Aldehydes 

Ketones 

Quinones 

Carboxylic  acids 
Acid  anhydrides 
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2  3  4  5  6  7  8 
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Figure  4.  Occurrence  frequency  of  the  principal  chealcal  groups  In  8462 
organic  coapounds  and  their  derivatives.  A  compound  aay  Include  functional 
groups  listed  above  It,  but  not  those  listed  below. 
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I.  TABULATION  OF  LASER-BREAKDOWN  THRESHOLDS  OF  SOLIDS 


Even  though  the  present  concern  with  low- loss  fibers  is  reducing  the 
scattering  and  absorption  in  order  to  increase  the  useful  length  of  fibers, 
laser  damage  of  fibers  is  of  interest  in  at  least  some  applications.  As 
seen  in  the  compilation,  only  a  small  fraction  of  the  laser-damage  studies 
were  specifically  concerned  with  fibers,  but  the  general  damage  results 
are  also  of  interest  in  the  study  of  fiber  of  course. 

There  have  been  so  many  experimental  investigations  of  laser-damage 

1  2 

of  solids  in  recent  years  *  that  a  tabulation  of  the  results  should  be 
useful.  Smith  tabulated  a  number  of  breakdown  thresholds  and  discussed 
the  theories  of  breakdown.  The  present  paper  is  a  tabulation  of  recent 
work  in  laser-induced  damage.  The  principle  source  of  information  for 
this  tabulation  was  the  annual  publication  by  the  U.S.  Department  of 
Commerce  and  the  National  Bureau  of  Standards  entitled  11  Laser- Induced 
Damage  in  Optical  Materials."^  The  purpose  of  the  tabulation  is  to 
provide  an  overview  of  the  vast  amount  of  data  in  the  literature  and 
to  identify  trends  in  the  data. 

The  complete  tabulation,  which  is  rather  long,  will  be  presented 
elsewhere.4  The  summary  figures  are  included  here  as  Figs.  1  through 
27.  The  following  materials  are  included  in  the  complete  tabulation,  but 
deleted  in  the  sunmary  figures  because  only  a  few  points  were  available: 

KF,  Rbtt,  NaBr,  Nal,  KI,  Rbl,  Csl,  BeF2,  CdTe,  GaAsP,  GaP,  ZnS,  Ba2NaNb015, 
Gd2Ga5°i2,  KH2P04,  La2Be205,  LUO3,  LiNb03,  LiTa03,  MgA£204,  MgO,  NH4H2P04 
(ADP),  SrTi03,  Ti02,  YAC03,  Y3AC5012,  Ag3AsS3,  As2S3,  BK-7  glass,  diamond. 


Sec.  I 


195 


CdCrgS^,  Cervit  glass,  KRS-5,  KRS-6,  LiYF4,  and  various  glasses. 

The  most  obvious  feature  of  the  results  is  large  spread  in  the 
experiment  results  for  the  same  material  at  the  same  wavelength.  Possible 
sources  of  the  great  variation  in  results  include  the  most  likely  source— 

•  extrinsic  sources  of  damage  that  are  different  in  different 
samples— 

and  the  following  other  sources: 

•  inaccuracies  in  the  irradtance  resulting  from  inaccuracies  in  the 
focal-spot  size,  particularly  that  resulting  from  aberrations  in 
the  optical  system, 

•  spatial-  and  temporal -spikes  (regions  of  high  trradiance)  in  the 
laser  beam  resulting  from  mode  locking  or  multimode  operation, 
for  example, 

•  focal-volume  dependence  of  the  breakdown  energy  density  resulting 
from  fundamental  effects  (an  unlikely  source), 

•  focal -volume  dependence  of  the  breakdown  energy  density  resulting 
from  technical  effects  (such  as  a  greater  probability  of  having 
an  easily  damaged  imperfection  in  the  focal  volume),  and 

•  effects  of  different  pulse  shapes,  both  as  functions  of  time  and 
space. 

The  effect  of  including  extrinsic  sources  of  damage,  aberration  in  the 
optical  system,  and  easily  damaged  imperfections  in  large  beams  is  to  reduce  the 
breakdown  energy  density;  whereas  the  effect  of  regions  of  high  irradiance 
(spikes)  is  to  increase  the  breakdown  energy  density.  It  is  likely  that  the 
larger  values  of  the  breakdown  energy  density  are  closer  to  the  correct  intrinsic 

I 
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values  because  the  samples  without  extrinsic  sources  of  damage  have 
greater  values  of  the  breakdown  energy  density.  Thus,  the  larger  values 
tend  to  be  more  representative  of  intrinsic  results.  This  tendancy  is  not 
universally  true  and  can  be  misleading  because  larger  experimantal  values 
can  also  be  ^caused  by  spiking  and  possibly  by  other  effects  that  have 
been  overlooked. 

Trends  in  the  data  are  difficult  to  identify  because  the  spread 

in  the  values  is  so  great.  Thus,  the  following  trends  should  not  be 

considered  as  fundamental  or  indicative  of  trends  expected  for  a  given 

controlled  experiment.  The  breakdown  energy  density  tends  to  increase 

with  increasing  pulse  duration.  The  spread  in  data  is  too  great  to 

determine  a  meaningful  relation  between  Itp  and  t  ,  but  the  relations 
1/2 

Itp  ~  tp  '  and  Itp  -  tp  are  sketched  for  copper  at  10.6  ym  in  Fig.  9, 
for  fused  silica  at  1.06  ym  in  Fig. 22,  and  for  sodium  chloride  at 

1/2 

10.6  ym  in  Fig.  19.  For  copper  and  fused  silica,  the  relation  Itp~  tp 
is  a  better  fit  to  the  data;  whereas,  for  sodium  chloride,  Itp  -  tp  is 
a  better  fit.  The  great  spread  in  the  data  prevents  a  convincing  choice 
between  even  the  tp^2  and  tp  trend. 

The  damage  energy  densities  at  2.7  ym  and  3,8  ym  tend  to  be  greater 
than  at  the  other  wavelengths.  However,  this  tendency  is  largely 
accounted  for  by  the  longer  pulse  durations  at  2.7  ym  and  3.8  ym.  Other¬ 
wise,  there  is  not  a  convincing  trend  of  Itp  as  a  function  of  wavelength 
on  the  region  from  0.53  ym  to  10.6  ym.  The  lack  of  a  striking  wavelength 
dependence  does  not  imply  that  the  intrinsic  breakdown  mechanism  is  wave¬ 
length  Independent  because  the  data  here  is  believed  to  be  largely  deter¬ 
mined  by  extrinsic  processes. 
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In  almost  all  of  the  pulsed-laser  damage  studied  to  date,  for  values 
of  irradiance  near  the  damage  threshold  the  damage  occurs  at  small  Isolated- 
spots  within  the  laser  beam  rather  than  uniformly  over  the  beam.  The 
known  possible  sources  of  the  isolated-spot  damage  are  all  extrinsic. 

Thus,  extrinsic  mechanisms  for  damage  and  for  plasma  initiation  dominate 
intrinsic  processes.  In  fundamental  experimental  investigations, 
intrinsic  breakdown  is  studied  by  focusing  a  laser  beam  into  a  small 
volume  of  the  bulk  material  away  from  the  surface  of  the  material.  Many 
measurements  are  made,  and  those  with  low  damage  thresholds  are  ignored. 
There  still  is  no  way  to  verify  that  the  large  thresholds  are  intrinsic. 

A  dramatic  illustration  of  this  uncertainty  is  the  new  high  damage- 

5 

threshold  data  from  Manenkov  and  coworkers  reported  at  the  1977 
Boulder  Damage  Symposium,  If  corraborated,  their  new  data  implies 
that  most  previous  damage  thresholds  measured  at  1.06  pro  were  extrinsic, 
rather  than  intrinsic  as  previously  believed. 

One  of  the  most  important  results  form  laser-damage  studies,  which 
if  often  overlooked,  is  that  surface-damage  thresholds  are  consistently 
low  when  large  beams  are  used.  For  small  beams,  it  is  possible  to 
obtain  larger  thresholds  for  some  shots— those  which  miss  the  low-damage- 
threshold  isolated  spots.  For  damage  by  large  laser  beams,  the  single¬ 
pulse  damage  thresholds  for  various  pulse  lengths  are  as  follows: 

2 

- — 1-10  0/cm  ,  for  -lOpsto  -10ns  pulses 

2 

few  hundred  0/cm  ,  for  microsecond  pulses  , 

These  values  are  the  same  for  metals,  dielectrics,  and  coatings. 


For  comparison,  the  Intrinsic,  electron-avalanche-breakdown 

thresholds  for  alkali  halides  at  1.0  ym  and  10ns  are  approximately 
3  2 

1.5  x  io  0/cm  ,  which  is  two- to  three-orders  of  magnitude  greater  than 
most  of  the  large-beam  surface  damage  values  above. 
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FIGURE  CAPTIONS  FOR  SEC.  I  j 

i 

t 

i 

Figs.  1-27.  Summaries  of  damage  thresholds  for  materials  with  more  than 
a  few  measurements.  Very  well  characterized  experiments 
are  denoted  •  or  x,  and  less  well  characterized  experiments 
are  denoted  o  or  a. 
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J.  OPTIMUM  FREQUENCY  RANGE 


In  this  section  it  is  shown  that  the  choice  of  the  optimum  frequency 

between  the  reststrahlen-  and  electronic -absorption-edges  involves  compromises. 

The  minimum  extrinsic  absorption  is  likely  to  be  near  1  pm.  The  zero- 

2  2 

dispersion  crossover,  as  the  frequency  at  which  d  n/dX  =  0  is  called, 
occurs  at  longer  wavelengths  for  most  materials  of  interest.  Also,  the 
scattering  is  likely  to  be  excessive  at  1  ym.  Operation  at  frequencies 
below  the  fundamental  reststrahlen-absorption  peak  should  be  considered. 

The  V-curves  of  Sec,  G  are  useful  in  establishing  the  range  of  wave¬ 
lengths  over  which  a  material  can  be  used  without  being  limited  by  intrinsic 
absorption  or  intrinsic  scattering.  The  choice  of  the  optimum  wavelength 
within  this  region  depends  on  such  technical  factors  as  reducing  the 
extrinsic  scattering  and  extrinsic  absorption. 

The  optimum  wavelength  region  between  the  reststrahlen-absorption 
peak  and  the  electronic-absorption  edge  for  minimum  extrinsic  absorption 
is  expected  to  be  near  1  pm.  The  absorption  peaks  of  molecular-ion 
impurities  occur  at  longer  wavelengths,  as  seen  in  Figs.  1-4  of  Sec.  H. 
Overtones  can  occur  in  the  1  pm  region,  but  overtone  absorption  is  weaker 
than  the  fundamental  absorption.  Electronic  absorption  is  strong  at 
wavelengths  shorter  than  1  pm.  See  the  tabulation  of  visible-  and  infrared- 
impurity-absorption  spectra  by  Sparks,  Vora  and  Flannery.1  The  lower 
absorption  at  1.06  pm  than  at  2.7,  3.8,  5.25,  and  10.6  pm  is  seen  in  the 
figures  in  Sec.  F. 

The  zero-dispersion  crossover,  in  fiber-optics  terminology,  is  the 

p  p 

wavelength  at  which  d  n/dX  -  0.  The  term  dispersion,  when  referring  to 


Sec.  J 


202 


the  index  of  refraction  n,  has  at  least  three  meanings— the  general  wave¬ 
length  dependence  of  n;  the  values  of  dn/dX;  and  the  value  of  d^n/dX2. 

In  order  to  avoid  confusion,  the  dispersion  of  fiber  optics,  which  is 
?  7 

d  n/dx  ,  will  be  called  material  dispersion.  As  explained  in 

Sec.  A,  the  material  dispersion,  along  with  the  inherent  dispersion  of  the 

mode  that  propagates  in  the  fiber,  determines  the  spatial  spreading  of  the 

different  frequency  components  of  the  information-carrying  laser  signal. 

2  2 

The  frequencies  at  which  d  n/dX  is  zero  are  indicated  on  the  V  curves 

in  Sec.  G.  The  values  were  obtained  from  the  minima  of  plots  of  dn/dX. 

Such  materials  as  fused  silica,  lithium  fluoride,  calcium  fluoride,  and 

sodium  fluoride  have  zero-dispersion  crossovers  between  1  and  1.5  pm. 

Covalently  bonded  materials  consisting  of  light  elements  tend  to  have 

their  zero-dispersion  crossovers  at  such  short  wavelengths,  while 

ionically  bonded  materials  consisting  of  heavy  elements  tend  to  have  their 

zero-dispersion  crossovers  at  longer  wavelengths. 

2  2 

For  the  materials  having  d  n/dX  *  0  at  short  wavelengths, the  extrinsic 
absorption  and  material  dispersion  can  both  be  small  in  this  wavelength 
region  near  1  pm.  Unfortunately,  the  scattering  tends  to  be  large  at  such 
shorter  wavelength,  at  which  the  extrinsic  absorption  and  dispersion  tend 
to  be  greater  for  these  materials.  For  such  other  materials  as  NaCl,  KC1, 

KBr,  KRS-5,  KRS-6,  AgCl ,  KI ,  CsBr  and  CSI  that  have  d2n/dX2  =  0  at  longer 
wavelengths  (between  3  pm  and  6  pm  for  the  materials  listed),  the  extrinsic 
absorption  is  likely  to  be  greater  at  the  longer  wavelength  of  zero-dispersion 
crossover. 

Simultaneously  achieving  low  material  dispersion,  low  extrinsic 
absorptance,  and  low  extrinsic  scattering  will  require  compromises  in  the 
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choice  of  materials  and  wavelengths.  The  intrinsic  absorption  and 
intrinsic  scattering  must  be  low  also,  but  intrinsic  scattering 
and  absorption  are  expected  to  be  negligible  with  respect  to 
extrinsic  scattering  and  absorption  in  most  materials  of  interest,  as 
discussed  in  Sec.  6. 

Consider  an  example.  Lithium  fluoride  has  zero  material  dispersion 
at  X  =  1.2  pm.  At  this  wavelength,  the  intrinsic  absorption  is  entirely 
negligible,  as  seen  in  the  figures  in  Secs.  F  and  G,  and  the  Brill ouin- 

_3 

scattering  coefficient  is  gg  =  2  x  10  dB/km.  Thus,  lithium  fluoride 
is  a  candidate  material,  based  on  zero-dispersion  crossover  and  intrinsic 
absorption  and  scattering,  even  though  the  zero-dispersion  crossover  at 
-1.2  pm  is  well  away  from  the  minimum  of  the  V  curves  at  X  =  2.1  pm 
(from  Sec.  G).  Such  fundamental  considerations  only  establish  that 
lithium  fluoride  is  a  reasonable  candidate.  The  more  difficult  task 
of  determining  whether  the  extrinsic  scattering  and  extrinsic  absorption 
can  be  reduced  to  tolerable  levels  involves  technical  factors,  rather 
than  fundamental  ones. 

another  example  of  material  candidates,  neither  TI  1173  glass  nor 

O 

gallium  arsenide  would  be  good  candidates  for  fibers  with  g_v  *  10'  dB/km 

C/\ 

_2 

because  intrinsic  absorption  and  scattering  process  give  g  >  10  dB/km 

vA 

_? 

for  TI  1173  glass  and  g  =  10  dB  at  the  minimum  of  the  V  curve  at 

CA 

X  =  6  pm  for  gallium  arsenide. 

Finally,  operation  at  frequencies  below  the  fundamental  reststrahlen 
absorption  peak  should  be  considered.  Extrinsic-  and  intrinsic-scattering 
and  extrinsic  absorption  should  be  much  less  severe.  It  must  be  determined 
whether  sufficiently  low  frequencies  to  avoid  excessive  intrinsic  absorption. 
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especially  at  room  temperature  and  above,  can  be  used.  Other  fundamental 
and  technical  limitations  should  be  carefully  considered. 
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K.  CHOICE  OF  GLASS-,  POLYCRYSTALLINE-, 
OR  SINGLE-CRYSTAL-FIBERS. 


Silica  fibers  are  believed  to  be  at  the  intrinsic -loss  limit,  after 
many  years  of  investigation.  New  materials  will  be  required  in  order  to 
attain  a  loss  less  than  the  lowest  silica-fiber  value  of  0.3  dB/km. 

It  is  too  early  in  the  low-loss-fiber  program  to  choose  between  glass-, 
polycrystalline-,  and  single-crystal-fibers.  The  central  issues  that  must 
be  settled  in  order  to  make  the  choice  are  identified  in  the  present  section. 

The  overwhelming  factor  at  present— the  ease  of  fabrication  of  glass 
fibers — favors  glass  fibers.  Any  new  glass-fiber  candidate  material  must 
be  fabri eatable  by  drawing  or  by  another  process  that  avoids  the  current 
difficulties  of  fabricating  crystalline  fibers,  or  the  advantages  of  the  glass 
surely  will  be  lost.  This  is  especially  important  because  some  candidate 
glasses  apparently  cannot  be  drawn.  Furthermore,  the  manufacturing  of 
significant  quantities  of  some  glasses  that  can  be  made  in  small  amounts 

may  be  impractical  by  current  techniques. 

If  the  ultimate  intrinsic-operation  limits  could  be  attained,  crystalline 
fibers  would  be  favored.  In  particular,  the  limited  number  of  glasses  , 
the  greater  intrinsic  scattering  in  glasses  compared  to  crystalline  materials, 
and  the  lower  intrinsic  absorption  of  many  crystals  favor  crystalline 
fibers.  Even  though  there  are  several  candidate  glasses,  the  small 
number  of  candidates  is  evinced  by  the  statement1  "...that  ZnU2 
glass,  though  water  soluble,  is  a  reasonable  and  possibly  unique  choice 
for  meeting  the  above  objectives  [of  attaining  ultralow  losses]."  Single¬ 
component  glasses  are  believed  to  be  required  to  avoid  excessive  scattering. 

In  order  to  reap  the  potential  benefits  of  crystalline  fibers,  the  strong 


4 


scattering  by  several  extrinsic  processes  must  be  greatly  reduced.  Scattering 
by  both  bulk-  and  surface-imperfections  in  current  crystalline  fibers,  both 
polycrystalline  and  single-crystal,  are  shown  in  Secs.  B  through  E  to  be 
extremely  strong.  A  new  fabrication  method  is  needed  for  crystalline  fibers. 
Removing  the  surface  and  bulk  scattering  should  be  the  main  goal  of  fabrication 
programs . 

The  problem  of  whether  single-crystal  fibers  will  be  required  or 
whether  polycrystalline  fibers  will  be  satisfactory  is  currently 
overshadowed  by  the  problem  of  whether  any  crystalline  fiber— single¬ 
crystal  or  polycrystalline— can  have  sufficiently  low  scattering. 

There  are  several  sources  of  loss  in  polycrystalline  fibers  that  do 
not  exist  in  crystalline  fibers.  It  is  well  known  that  contaminated  grain 
boundaries  are  a  major  source  of  absorption  loss  in  bulk  polycrystalline 
materials.  The  original  Kodak  Irtran  materials  afford  good  examples  of 
high  loss  resulting  from  contaminated  grain  boundaries.  However,  it  was 
established  in  high-power  laser-window  studies  that  internally  formed  grain 

boundaries  can  be  sufficiently  clean  to  have  negligible  absorption,  at  least 

-4  -1 

at  the  level  of  Bab  =10  cm  .  Thus,  grain- boundary  contamination  is  not 

necessarily  a  problem,  but  it  is  not  known  whether  polycrystalline  fibers 

can  be  manufactured  with  sufficiently  clean  grain  boundaries  to  avoid 

-8  -1 

excessive  absorption  and  scattering,  especially  at  ggx  ?  10  cm"  . 

An  important  method  of  strengthening  materials  is  to  use  polycrystalline 
materials  and  reduce  the  grain  size.  Thus,  small-grain  polycrystalline  fibers 
are  expected  to  be  stronger  than  single-crystal  fibers.  Grain  growth  is 
likely  to  be  a  problem.  However,  alkali-halide  crystals  grown  by  the  reactive- 
atmosphere  process  are  free  of  OH "ions  and  do  not  show  grain  growth.  The 
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strength  of  polycrystalline-  and  single-crystal  fibers  was  considered  by 

2 

Harrington,  Turk,  Henderson,  and  ttyer, 

3 

Fredricks  pointed  out  that  single-crystal  alkali-halide  fibers  may 
have  a  characteristic  that  will  limit  their  usefulness.  When  alkali  halide 
crystals  bend,  they  do  so  by  sliding  on  the  <100>  planes,  thereby  forming 
a  set  of  steps  around  the  corner.  Also,  once  bent,  the  fiber  is  likely 
to  break  on  the  next  attempt  to  bend  it.  The  bonding  between  adjacent 
<100>  planes  is  markedly  weaker  than  any  other  planes  in  the  sodium-chloride 
structure.  The  steps  formed  in  typical  bending  are  large  with  respect  to 
the  wavelength  of  the  light.  The  resulting  scattering  will  be  intolerably 
great,  as  shown  in  Sec.  B.  Fittings  could  be  used,  with  straight  sections 
of  fiber  in  between  fittings,  but  this  solution  would  surely  be  impractical 
in  many  applications.  Installation  would  be  a  problem.  The  bending 
effect  would  limit  the  size  of  the  take-up  reels  used  in  producing  the 
fibers. 

Processes  other  than  extrusion  are  of  current  interest  because 

current  extrusion  processes  have  been  unsuccessful,  Fredricks  has  made 

suggestions  to  improve  extruded  fibers.  First,  it  may  be  possible  to 

extrude  crystalline  fibers  with  smooth  surfaces.  Surface  layers  of 

2 

current  extruded  potassium-chloride  fibers  split  from  the  main  fiber 
in  a  rather  regular  pattern,  somewhat  like  fish  scales  or  like  a 
banana  whose  skin  had  been  partially  peeled  so  that  the  peels  form  an 
acute  angle  with  the  meat.  Fibers  are  made  by  extruding  a  potassium- 
chloride  crystal  through  a  parabolic  orifice  diamond  die  at  300-400°C. 

The  potassium-chloride  fiber  must  have  <100  planes  parallel  to  Its 
axis  and,  of  course,  transverse  to  it.  This  puts  the  weak  crystallographic 


directions  in  the  proper  orientation  to  form  the  surface  scales. 

If  the  surface  tensions,  y,  of  potassium-chloride,  diamond, and  the 
surrounding  atmosphere  are  such  that  the  Inequality 

y  (diamond  -  gas)  +  y  (diamond  -  KC1)  <  y  (KC1  -  gas). 

Is  valid,  the  potassium-chloride  will  wet  the  diamond.  In  addition,  a 
second  attractive  force  between  potassium-chloride  and  the  diamond  die  can 
generatedlf  the  hydrostatic  pressure  used  to  extrude  the  fiber  is  sufficient 
to  reduce  the  <100>  steps,  which  separate  the  potassium-chloride  -  diamond 
surfaces, to  less  than  100  nm.  This  Is  an  induced-dipole  force, or  van  der 
Waals  force.  A  common  example  is  the  attractive  force  between  two  Johnson 
blocks.  These  highly  polished  blocks  cannot  be  easily  pulled  apart,  but 
must  be  slid  along  the  polished  surfaces  to  be  separated. 

Data  on  surface-tension  measurements  between  diamond  -  potassium-chloride 
gas  systems  has  not  yet  been  found.  Even  If  they  were  available  they  would 
have  to  be  used  with  caution  as  trace  water  changes  the  properties  greatly. 

The  Induced -dipole  forces  will  be  a  function  of  the  extrusion  force  and 
roughness  of  the  die. 

As  the  potassium-chloride  crystal  Is  extruded  through  the  die,  it 
becomes  strongly  bound  to  the  diamond.  As  It  exits  the  parabolic  orifice 
It  tends  to  follow  the  surface  of  the  die  until  the  fiber  splits  on  the 
<100>  cleavage  planes  parallel  to  the  fiber  axis  and  then  breaks  on  the 
transverse  <100 >  cleavage  planes.  Since  both  surface  tension  and  induced- 
dipole  forces  tend  to  retard  the  surface  planes  when  the  surface  “peels" 
back,  the  scales  bend  slightly  away  from  the  fiber  as  they  attempt  to 
relieve  the  Internal  strain.  There  must  be  a  minimum  of  four  peels  around 
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the  fiber.  More  can  be  present  If  grain  boundaries  occur  across  the  fiber 
axis. 

This  scaling  effect  can  be  avoided  by  reducing  the  attractive  forces 
between  potassium-chloride  and  the  die.  The  surface-tension  problem  can 
be  avoided  by  using  dry  potassium-chloride  and  a  quartz  die.  The 
van  der  Waals  attraction  can  be  reduced  by  lowering  extrusion  pressures, 
and  the  tendency  for  splitting  lessened  by  removing  the  parabolic  exit  on 
the  die.  The  exit  should  be  cut  square  or  at  a  small  angle  above  the  exit 
hole  of  the  die. 

Molten  potassium  chloride  could  be  extruded  through  a  quartz  orifice 
and  either  recrystallized  as  one  does  for  wire,  or  for  a  more  practical 
fiber  quenched  to  a  glassy  state.  A  rather  simple  apparatus  could  be 
constructed  from  a  quartz  crucible  with  a  quartz  capillary  attached  to  Its 
bottom.  The  capillary  Is  drawn  to  fiber  size  and  cut  square.  The  top  of 
the  crucible  and  the  tube  enclosing  the  orifice  must  be  sealed  to  allow 
hydrogen  chloride  flushing.  After  flushing, the  potassium  chloride  would 
be  quenched  by  cooled  argon  blown  against  the  orifice  and  down  the  enclosing 
tube  around  the  fiber.  If  these  fibers  have  satisfactorily  optical 
properties,  a  take-up  spool  could  be  added  to  obtain  longer  samples. 

There  is  a  possibility  that  Fredricks  could  carry  out  an  experimental 
program  to  obtain  Improved  extrusion  at  Oregon  State  University  If  there 
Is  Interest. 
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M.  SUMMARY  &  CONCLUSIONS 


After  many  years  of  research,  fused-silica  fibers  are  believed  to  be 
at  their  intrinsic-loss  limit,  which  is  an  extinction  coefficient  of 
Bgx  =  0.3  dB/km  at  the  minimum  of  the  V  curves  at  a  wavelength  near  1.6  pm. 
Thus,  new  materials  will  be  required  in  order  to  attain  a  lower  loss. 

The  V  curves,  as  the  plots  of  the  intrinsic-scattering  coefficient 
and  intrinsic-absorption  coefficient  are  called,  are  useful  in  coarse 
screening  of  candidate  fiber  materials  and  in  establishing  the  useful 
wavelength  region  of  a  given  material,  as  illustrated  in  Fig.  8  of 
Sec.  G.  For  low-loss  fibers,  with  extinction  coefficients  fi  «  0.1  dB/km 

v.  i \ 

=  2.3  x  10~7  cm”\  both  the  scattering  and  absorption  for  wavelengths 
within  this  useful  wavelength  region  are  determined  by  extrinsic  processes 
in  most  materials  of  interest. 

In  communications  applications,  it  is  usually  necessary  to  operate 

within  a  narrow  wavelength  region  near  the  zero-di spersion-crossover 

2  2 

wavelength  which  is  the  wavelength  at  which  dn  /d\  =  0,  where  n 

is  the  real  part  of  the  index  of  refraction.  Thus,  the  problem  of 
obtaining  a  satisfactory  fiber  material  is  as  follows:  (1)  Find  a 
material  with  X^  near  the  operating  wavelength  (or  conversely,  choose 
the  operating  wavelength  equal  to  Xq^  for  the  most  promising  candidate 
material,  if  possible).  (2)  Verify  that  the  intrinsic  loss,  given  on 
the  V  curves  in  Sec.  G,  is  at  or  below  the  required  value.  It  is  not 
necessary  that  the  minimum  of  the  V  curve  fall  near  Xq^.  It  is  sufficient 
that  is  in  the  useful  wavelength  region.  (3)  Improve  the 
material  purity  and  fiber-fabrication  process  to  obtain  the  required 
low  value  of  absorption  and  scattering. 
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This  last  step  of  reducing  the  extrinsic  absorption  and  extrinsic 
scattering  is  expected  to  be  difficult  for  low-loss  fibers.  In  Secs.  B 
through  E,  it  is  shown  that  any  one  of  several  scattering  mechanisms  can 
give  rise  to  a  scattering  coefficient  that  is  five  to  eight  orders  of 
magnitude  greater  than  the  typical  goal  of  3  =  10  dB/km.  In  Sec.  H 

Ca 

-4 

it  is  shown  that  impurities  at  the  level  of  approximately  10  parts  per 

billion,  which  is  too  low  to  even  measure  for  important  impurities,  can 

_2 

give  rise  to  an  absorption  coefficient  of  3ab  =  10  dB/km. 

The  success  of  the  low-loss  fiber  program  depends  on  the  ability  to 
solve  the  technical  problems  of  purifying  materials  and  maintaining  their 
purities  during  fabrication,  installation,  and  field  use  in  order  to 
obtain  low  extrinsic  absorption  and  scattering  and  of  fabricating  fibers 
with  low  scattering  loss  that  will  remain  low  during  installation  and 
operation.  These  severe  problems  are  in  addition  to  usual  problems  of 
obtaining  claddings,  achieving  index  gradation,  coupling  fiber  sections 
together,  and  obtaining  sufficiently  strong  fibers.  This  pessimistic 
picture  establishes  the  priorities  in  the  low-loss  fiber-optics  program, 
rather  than  suggesting  that  the  program  cannot  succeed.  After  all,  a 
loss  of  ~0.3  dB/km  has  been  attained  in  spite  of  the  difficulties. 

It  is  too  soon  in  the  low-loss  fiber  program  to  choose  between  glasses, 
polycrystals,  and  single  crystals.  The  problem  of  choosing  between  single¬ 
crystal  fibers  and  polycrystalline  fibers  is  overwhelmed  at  present  by  the 
problem  of  whether  any  crystalline  fiber  can  have  sufficiently  low 
scattering.  The  choice  between  polycrystalline-  and  single-crystal-fibers 
cannot  be  made  until  the  excessive  scattering  problem  of  crystalline  fibers 
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is  solved.  Nevertheless,  the  issues  involved  in  the  choice  between  single¬ 
crystals  and  polycrystals,  as  well  as  the  choice  between  glasses  and 
crystals,  are  discussed  in  Sec.  L. 

Ease  of  fabrication  is  an  overwhelming  factor  in  favor  of  glasses 
that  can  be  drawn  as  successfully  as  fused  silica.  The  limited  number 
of  glasses,  the  greater  intrinsic  scattering  and  absorption  in  glasses  than 
in  crystals,  and  the  difficulty  in  manufacturing  large  Quantities  of  some 
glasses  favor  crystalline  fibers.  It  ic  possible  that  single-crystal 
alkali -halide  crystals  may  not  be  useful,  even  if  satisfactorily  fabri¬ 
cated,  because  bending  the  fiber  would  result  in  surface  imperfections 
(<‘i00>  "steps")  that  would  cause  excessive  scattering,  and  multiple 
bendings  may  cause  fracture.  An  improved  extrusion  method  suggested 
by  Fredricks  and  discussed  in  Sec.  L  possibly  could  result  in  greatly 
improved  extruded  fibers. 

The  choice  of  the  optimum  frequency,  for  frequencies  between  the 
reststrahlen-  and  electronic-absorption  edges,  involves  compromises. 

The  minimum  extrinsic  absorption  is  likely  to  be  near  1  pm.  The  zero- 
dispersion  crossover  occurs  at  wavelengths  longer  than  1  pm  for  most 
materials  of  interest.  Also,  the  scattering  is  likely  to  be  excessive 
at  1  pm.  Operation  at  frequencies  below  the  fundamental  reststrahlen- 
absorption  peak  should  be  considered. 

Intrinsic-Brillouin  and  intrinsic-P.ayleigh  scattering,  which 
determine  the  lowest  possible  scattering  that  can  be  attained  in  fiber 
materials,  have  typical  scattering  coefficients  B$c  =  Bg  +  ^ 

of  1  x  10-9  cm’1  =  4  x  10'4  dB/km  for  lithium  flouride  to  8  x  io~7  cm’1  = 

0.4  dB/km  for  zinc  selenide,  both  at  wavelength  A  =  1  pm.  The  total  intrinsic 
scattering  coefficient,  for  Brillouin  and  intrinsic  Raman  scattering. 


can  be  much  greater,  by  a  factor  of  40  in  the  examples  considered, than 
the  Brillouin  scattering,  in  contrast  to  the  tacit  assumption  made  in  the 
the  standard  practice  of  using  the  Brillouin  scattering  as  that  of  the 
V  curves  for  crystalline  materials.  Electronic  scattering  from 
vacancies  or  impurity  ions  that  are  either  isolated  or  in  small 
clusters  is  negligible,  but  electronic  scattering  from  extended 
imperfections,  such  as  precipitates,  grain  boundaries,  Suzuki  phases, 
and  dislocations  can  be  great. 

The  substantial  variations  in  both  the  laser-breakdown  thresholds 
and  the  amount  of  scattering  that  is  observed  as  various  regions  of 
crystal  are  probed  are  explained  in  terms  of  the  distribution  of 
impurities  and  imperfections,  which  change  in  time.  Quenching  temporarily 
reduces  clustering,  thereby  reducing  scattering  and  possibly  isolated- 
spot  laser  damage;  whereas  annealing  temporarily  increases  clustering, 
thereby  increasing  scattering  and  possibly  isolated-spot  laser  damage. 
Clusters  are  expected  to  form  in  time  even  in  quenched  crystals,  and 
0H~  and  H^O  substantially  penetrate  alkali-halide  crystals  if  they 
contain  a  few  parts  per  million  of  divalent  ions. 

The  substantial  variations  in  both  the  laser- breakdown  thresholds 
and  the  amount  of  scattering  that  is  observed  as  various  regions  of 
crystal  are  probed  can  be  understood  as  follows:  Even  though  the 
brief  review  in  Sec.  C-III.  of  the  current  state  of  impurity 
chemistry  in  the  simplest  ionic  solids,  the  alkali  halides,  is  far 
from  extensive,  it  should  disspell  the  common  view  that  simple  point 
defects  are  the  predominate  impurity  species  in  alkali  halides.  Further¬ 
more,  most  textbooks  create  the  impression  that  these  "simple  point 


defects"  are  frozen  in  some  reasonably  homogeneous  distribution  at  room 
temperature..  There  is  ample  evidence  that  this  is  not  true. 

In  Sec.  D  it  is  shown  that  the  scattering  by  the  strain  field  around 
an  impurity  ion  is  negligible  and  cannot  account  for  the  well-known 
difference  between  the  ion  radius  and  the  effective  radius  for  scattering 
In  Sec.  E  it  is  shown  by  developing  a  multiple-scattering  formalism  and 
applying  the  results  to  scattering  by  randomly  located  spherical  voids 
that  the  multiple-scattering  reduction  of  the  individual  scattering 
results  is  entirely  negligible  in  low-loss  fibers  because  the  scattering 
must  be  small  and  multiple-scattering  effects  are  negligible  for  small 
scattering. 

The  lowest  experimental  values  of  infrared  optical -absorption 

coefficients  from  the  literature  are  tabulated  in  Sec.  F  as  plots  of 

absorption  coefficients  as  functions  of  frequency.  Both  intrinsic  and 

extrinsic  values  are  tabulated.  The  intrinsic,  multiphonon  absorption 

curves  must  be  extrapolated  in  order  to  obtain  the  low  values  of 

current  interest  because  the  state  of  the  art  in  reducing  the  absorption 

coefficient  of  window  materials,  and  even  in  measuring  the  absorption 
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coefficients,  limits  the  values  of  B  to  approximately  10  cm  to 
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10  cm  (44  to  4.4  dB/km).  These  extrapolated  curves  are  given  as  the 
absorption  parts  of  the  V  curves  in  Sec.  G.  Laser-damage  thresholds  in 
solids,  including  fibers,  are  tabulated  in  Sec.  I.  Publications  produced 
in  the  study  are  listed  in  Sec.  L. 
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